LNF-92/xxx 
CERN-TH.6504/92 
BUTP-92/20 



(N 

o\ 

^ SEMILEPTONIC KAON DECAYS 

4j5 : IN CHIRAL PERTURBATION THEORY 

J. Bijnens and G. Eckei^ 

o : CERN 

<N ; CH- 1211 Geneva 23 

00 

o 

On 



and 



P-i. 

i \ J. Gasser* 
Ph. 

<D , 

c-| . Institut fiir theoretische Physik 

Universitat Bern, Sidlerstrasse 5 

^ ■ CH-3012 Bern 



" Permanent address: Institut fiir theoretische Physik, Universitat Wien, Boltz- 

manngasse 5, A-1090 Wien, Austria 
* Research supported in part by Schweizerischer Nationalfonds 



CERN-TH.6504/92 
May 1992 



Abstract 



We present the matrix elements for the semileptonic kaon decays K^, K121+1-, K12V+V-, 
Ki 3 , Ki 3l and at next-to- leading order in chiral perturbation theory and compare the 
predictions with experimental data. Monte Carlo event generators are used to calculate the 
corresponding rates at DAFNE. We discuss the possibilities to improve our knowledge of the 
low-energy structure of the Standard Model at this and similar machines. 
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Note: 

The number of events quoted for DAFNE are based on a luminosity of 5 • 10 32 cm" 2 s _1 , 
which is equivalent [] to an annual rate of 9 • 10 9 (1.1 • 10 9 ) tagged (Kl) (1 year = 
10 7 s assumed). 

Whenever we quote a branching ratio for a semileptonic K° decay, it stands for the 
branching ratio of the corresponding Kl decay, e.g., 

BR(K° -> ir~l + v) = BR(K L -> tt ± / t z/) . 
More notation is provided in appendix [A|. 



1 P. Franzini, private communication. 
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Chapter 1 

INTRODUCTION TO CHIRAL 
SYMMETRY 



Chiral perturbation theory (CHPT) is a systematic approach to formulate the standard 
model as a quantum field theory at the hadronic level. In its general form, it uses only the 
symmetries of the standard model, in particular its spontaneously broken chiral symmetry. It 
is characterized by an effective chiral Lagrangian in terms of pseudoscalar meson fields (and 
possibly other low-lying hadronic states) giving rise to a systematic low-energy expansion of 
amplitudes [1], [|. 

In the formulation of Ref.0, one considers the generating functional Z[v, a, s,p] of con- 
nected Green functions of quark currents associated with the fundamental Lagrangian 

£ = C QCD + ?7 M (^ + 75<v)g - q(s - ilsp)q. (1) 

Cqcd i s the QCD Lagrangian with the masses of the three light quarks set to zero. The 
external fields i> M , a M , s and p are hermitian 3x3 matrices in flavour space. To describe elec- 
tromagnetic and semileptonic interactions, the relevant external gauge fields of the standard 
model are Q 

r tl = v ll + a ft = -eQA^ (2) 

/ M = ^-« M = —eQAn — 6 (W+T+ + h.c.) 

V2sin0v 

Q = idiag(2, -1, -1), T + = 

where the are Kobayashi-Maskawa matrix elements. The quark mass matrix 

M = di&g(m u ,m d ,m s ) (3) 

is contained in the scalar field s(x). The Lagrangian (j^) exhibits a local SU(3)l x SU(3)r 
symmetry 

1. 1. 
Q 9R^V- + lb)q + 9L^{l-75>q 

r» -> gRr^R + igRd^R 

lii -> gL^gl + WLd^gl (4) 
s + ip -> flr fl (a + ip)</J, 
g-i^i G SU(3) RjL . 

The generating functional Z admits an expansion in powers of external momenta and 
quark masses (CHPT). In the meson sector at leading order in CHPT, it is given by the 
classical action 

Z = J d 4 x£2(U, v , a, s,p). (5) 
1 We adopt the present conventions of the Particle Data Group gj. 
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C-2 is the non-linear a model Lagrangian coupled to the external fields v, a, s,p 

F 2 



where 



D ll U = d fl U-ir ll U + iUl ll , x = 277 (s + ip), 
and (A) stands for the trace of the matrix A. U is a unitary 3x3 matrix 

17*17=1, detC/ = l, 



(6) 
(7) 



which transforms as 



under SU(3)l x SU(3)r. U incorporates the fields of the eight pseudoscalar Goldstone 
bosons. A convenient parametrization is R 



[/ = exp (iV2$/F), $ 



??8 

\/2 Vq 



-7T 



V 



7T 



_ + 



(9) 



The parameters 7 1 and Bq are the only free constants at 0(p 2 ) 
constant in the chiral limit, 

= F(l + 0{m quark )) = 93.2MeV, 
whereas B is related to the quark condensate, 

(0|mt|0) = -F 2 77 (l + 0{m quark )). 



-K+ ^ 
-K° 

F is the pion decay 
(10) 

fin 



B always appears multiplied by quark masses. At 0(p 2 ), the product B m q can be expressed 
in terms of meson masses, e.g. 



Mi 



B (m u + m d ). 



(12) 



The Lagrangian @ is referred to as the effective chiral Lagrangian of 0(p 2 ). The chiral 
counting rules are the following: the field U is of 0(p°), the derivative <9 M and the external 
gauge fields t> M , are terms of 0(p), and the fields s,p count as 0{p 2 ). 

At order p 4 the generating functional consists of three terms : 

i) The one-loop graphs generated by the lowest order Lagrangian (Q). 

ii) An explicit local action of order p 4 . 



2 We follow the Condon-Shortley-de Swart phase conventions. 
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iii) A contribution to account for the chiral anomaly. 



We briefly discuss the contributions ii) and iii) and start with the local action of 0(p 4 ). 
It is generated by the Lagrangian £ 4 0: 

£ 4 = L^D^D^U) 2 + L 2 (D P U ] D V U){D»U ] D V U) 

+L 3 (D^D»UD u U i D»U) + U{D^D»U){^U + xU ] ) 
+L s (D lt rtD"UtfU + U^x)) + L,{ X ] U + X U ] ) 2 + L 7 &U - X U ] ) 2 
+L 8 (x f Ux^U + xtfxrf) - tL 9 (F^D,UD u U f + F£ U D^D U U) 
+L W (U^F^UF L ^) + L u (Vff + F LfJiV F£ u ) + L l2 { X ] x), (13) 

where 

Fg = d»r v -0V-i[rV] (14) 
F% v = dT-<9T-z[P,r]. 

The twelve new low-energy couplings Li, . . . , L\ 2 arising here are in general divergent (except 
L3, L7). They absorb the divergences of the one- loop graphs via the renormalization 

A = (4 7 r)-V- 4 {^ I -^(ln(47r) + r(l) + l)} (15) 

in the dimensional regularization scheme. The coefficients Tj are displayed in table |l|. They 
govern the scale dependence of the renormalized, finite couplings L^(fj), 

m»2)=m»i)+^h^- ■ (is) 

Observable quantities are independent of the scale /1, once the loop contributions are in- 
cluded. 

The constants F,B , together with L\, ...,L\ Q , completely determine the low-energy 
behaviour of pseudoscalar meson interactions to 0(p A ). L\ x and L\ 2 are contact terms 
which are not directly accessible to experiment. Similarly to F and B discussed above, the 
constants L\ are not determined by chiral symmetry — they are fixed by the dynamics of 
the underlying theory through the renormalization group invariant scale A and by the heavy 
quark masses m c ,mb, .... With present techniques, it is, however, not possible to evaluate 
them directly from the QCD Lagrangian. In the absence of such a calculational scheme, 
they have been determined by comparison with experimental low-energy information and by 
using large— arguments. The result is shown in column 2 of table [l], where L\, . . . , L r 1Q 
are displayed at the scale \x = M p . The experimental information underlying these values is 
shown in column 3. [Li, L 2 and L 3 are taken from a recent overall fit to K ei and tttt data ||, 
see also the section on decays in chapter g and Ref. M. L 4 , . . . , L 10 are from ||. For L 9 
see also |J. In Refs. |2|, [7| it was shown that the values for the L^(M p ) can be understood 
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Table 1: Phenomenological values and source for the renormalized coupling constants 
Ll(Mp). The quantities Tj in the fourth column determine the scale dependence of the 
LKn) according to Eq. (|i~6|). L r n and LT 12 are not directly accessible to experiment. 



i 


Ll(Mp) x 10 3 


source 


Ti 


1 


0.7 ± 0.5 


K e 4, TXTX — > TXTX 


3/32 


2 


1.2 ± 0.4 


K e 4, TXTX — > TXTX 


3/16 


3 


-3.6 ± 1.3 


K e 4, TXTX — > TXTX 





4 


-0.3 ± 0.5 


Zweig rule 


1/8 


5 


1.4 ± 0.5 


Fr '■ F-jx 


3/8 


6 


-0.2 ± 0.3 


Zweig rule 


11/144 


7 


-0.4 ± 0.2 


Gell-Mann-Okubo,L 5 , L 8 





8 


0.9 ± 0.3 


M K o - M K+ ,L 5 , 


5/48 






(2m s -m u - m d ) : (m d - m u ) 




9 


6.9 ± 0.7 


< r 2 > 7r 

em 


1/4 


10 


-5.5 ± 0.7 


TX — > eZ/7 


-1/4 


11 
12 






-1/8 
5/24 



in terms of meson resonance exchange. For recent attempts to evaluate Li directly from the 
QCD Lagrangian see ||.] 

Here, it is of interest to know which of the low-energy couplings occur in the matrix 
elements for the semileptonic kaon decays discussed in chapter [| This information is given 
in table 0. (There is an ambiguity concerning the bookkeeping of L 4 and L 5 : some of 
these contributions may be absorbed into the physical decay constants F n , F K . Here we 
have chosen the convention which corresponds to the amplitudes displayed in chapter [|. 
Furthermore, in decays, additional constants may occur via the form factor R which 
has not yet been worked out at one-loop level 0. This channel is therefore omitted in the 
table.) 

We now turn to point iii) above. A functional Z[U,l,r] which reproduces the chiral 
anomaly was first constructed by Wess and Zumino flTUfl . For practical purposes, it is useful 
to write it in the explicit form given by Witten fLT |: 

Z[U,l,r] wzw = -^/ M5 ^' Hm (SfS^^fS^) (17) 

W(U,l,r)^ = (Ul,lJ a U^rp + -Ul^rM a U^rp + iUd ll lJ a U^ p 
+id^rMlaU j rp - iZXtfraUlp + E L U^d v r a Ul p 
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Table 2: Occurrence of the low-energy coupling constants L\, . . . , L\o and of the anomaly in 
the semileptonic decays discussed in chapter ||. 







Kim 


K l3 


Kl3-y 


K+ -> 


K+ -> 

7T 7T e + Z/ e 


K° — > 

7T 7T _ e + Z/ e 












X 


X 














X 


X 














X 


X 


X 


U 










X 


X 




u 










X 


X 


X 


u 




X 


X 


X 


X 


X 


X 


Lg + Lio 


X 


X 




X 








Anomaly 


X 


X 




X 


X 




X 



-tfEtrfrJJlf, + tflMp + tfdjjp (18) 
—iEplJalp + -H^lvEJp - iE^E^E^) 

sj = t^c/ e* = era,//* 

A^c = 3 £oi23 = 1 

where (L <-> i?) stands for the interchange 

~ ^, ^ ~ r„, ^ «-> Ej. 

The integration in the first term in Eq. ([17]) is over a five- dimensional manifold whose 
boundary is four-dimensional Minkowski space, such that 

/ d 5 xe^ klm d m T ijkl = [ d 4 xe^T, vpa (19) 

JM 5 J 

according to Stoke's theorem. [This term involves at least five pseudoscalar fields and will 
not be needed in the following chapter.] The convention used in Eq. (|17|) ensures that 
Z[U, I, r] WZ w conserves parity and reproduces the anomaly under SU (3) l x SU(3)r transfor- 



mations in Bardeen's form |12| (in particular, it is invariant under transformations generated 
by the vector currents). 

The Wess-Zumino-Witten functional contains all the anomalies which contribute to the 
semileptonic meson decays considered in the following chapter. The relevant piece for e.g. 
decays is 

Z[U, I, r] wzw = £^ J d 4 x W < ^ffWSv* >+■■■ . (20) 



8 



This short introduction to CHPT (see Refs.il3|, EHl for more extensive treatments with 



references to the original literature) contains all the ingredients necessary for the calculation 
of semileptonic K decay amplitudes to 0(p 4 ) presented in the next chapter. For the low 
energies involved in these decays, the momentum dependence of the W propagator connecting 
to the lepton-neutrino pair in the final state can be neglected. The chiral realization of the 
non-leptonic weak interactions is discussed in the corresponding chapters on non-leptonic K 
decays. 
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Chapter 2 

SEMILEPTONIC KAON DECAYS 



1 Radiative K\i decays 



We consider the K[ 2l decay 

#+(p) - l + (pi)vi{pMq) [K l27 ] (1.1) 

where / stands for e or /i, and 7 is a real photon with q 2 = 0. Processes where the (virtual) 
photon converts into a e + e~ or /i + yU~ pair are considered in the next section. The K~ mode 
is obtained from ( |1 . 1| ) by charge conjugation. 



1.1 Matrix elements and kinematics 

The matrix element for K + — > / + z^7 has the structure 



T = -iG F eV:A {F K L» - H^l v } (1.2) 



with 



\2pq 2piq J 
T = u(p v )'f(l-<y B )v(p l ) 

= iV(W 2 )e^ aP q a p p - A(W 2 )(qWg^ v -WY) 
W» = {p-qT = (pi+ Pv y. (1.3) 

Here, e M denotes the polarization vector of the photon with g^e M = 0, whereas A, V stand 
for two Lorentz invariant amplitudes which occur in the general decomposition of the tensors 

= J dxe^ +iW y < I TV? m {x) IUM I K + (p) > , I = V, A . (1.4) 



The form factor A (V) is related to the matrix element of the axial (vector) current in (fE 
In appendix |T] we display the general decomposition of A^, V^ v for q 2 7^ and provide also 
the link with the notation used by the PDG [0 and in @, §]. 

The term proportional to in (|1.2|) does not contain unknown quantities - it is deter- 



mined by the amplitude of the nonradiative decay K + — > This part of the amplitude is 
usually referred to as "inner Bremsstrahlung (IB) contribution", whereas the term propor- 
tional to is called "structure dependent (SD) part" . 

The form factors are analytic functions in the complex H /2 -plane cut along the positive 
real axis. The cut starts at W 2 = (M K + 2M W ) 2 for A (at W 2 = (M K + M n ) 2 for V). In our 
phase convention, A and V are real in the physical region of decays, 

mf < W 2 < M 2 K . (1.5) 

The kinematics of (spin averaged) K\ 2l decays needs two variables, for which we choose 
the conventional quantities 

x = 2pq/M 2 K , y = 2p Pl /M 2 K . (1.6) 
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In the K rest frame, the variable x (y) is proportional to the photon (charged lepton) energy, 

x = 2E 1 /M K , y = 2E l /M K , (1.7) 
and the angle 6^ between the photon and the charged lepton is related to x and y by 

(1 - y/2 + A/2)(l - y/2 - A/2) = 

l-y/2 + A/2cos6 h yy K ' 

In terms of these quantities, one has 

W 2 = M 2 K {l-x) ; (q 2 = 0) . (1.9) 

We write the physical region for x and y as 

2^n< y < 1+r, 



or, equivalent ly, as 



l-\(y + A) <x< l-l( y -A) (1.10) 



< x < 1 — r\ 

\-x + -p—< y <l + r t (1.11) 
(1 -x) 



where 



1.1 • 1(T 6 (Z = e) 
4.6 • 1(T 2 (/ = /i) . 



1.2 Decay rates 

The partial decay rate is 



spins 

The Dalitz plot density 



dT = 2M^f £ m 2 <Wp;^,g). (1-13) 



dxdy 256vr 3 s ^ 



spins 



is a Lorentz invariant function which contains V and A in the following form 

= Pm(x,y) +p SD (x,y) + p mT (x,y) 

p m (x,y) = A m f m (x,y) 

Pso (x,y) = A SD M 2 K [(V + A) 2 f SD+ (x,y) + (V - A) 2 f SD -(x,y) 

y) = A mT M K [(V + A)f mT+ (x,y) + (V -A)f mT -(x,y)} (1.15) 
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where 



fm(x,y) 

fm- (?, y) 



l-y + ri 



x 2 + 2(1 - ar)(l - n) 



/lNT+ ' 
/lNT~ ' 



x,y) 
x,y) 



x 2 (x + y — 1 — r\ 
[x + y - 1 - ri] [(x + y - 1)(1 - x) - n\ 
[l-y + n] [(l-x)(l-3/)+r,] 
1 - y + n 



2xn{\ - n) 

x + y — 



a;(x + y — 1 — rj) 

l-y + ri 
x(x + y — 1 — ri) 



[(1 -x)(l -x-y)+r l ] 

(1 - x)(l - x -y) -r x 



x 



;i.i6) 



and 



A m — 

^SD — 
-^INT 

1 1 )! later convenience, we note that 

a 1 



4n 



A 



M, 



A 



-4, 



4r, 



32tt 2 
^a 



-M 



K 



M, 



A 



-4, 



.4. 



*7rr,(l-n) 2 VF K 



M, 



A" 



;i.i7) 



;i.i8) 



The indices IB, SD and INT stand respectively for the contribution from inner Bremsstrahlung, 
from the structure dependent part and from the interference term between the IB and the 
SD part in the amplitude. 

To get a feeling for the magnitude of the various contributions IBjSD 111 and INT 1 * 1 to the 
decay rate, we consider the integrated rates 



r 7 = / dxdypAx.y) ; I = SD ± , INT ± , IB 



1.19) 



where psd = Psd+ + Psd- etc. For the region Rj we take the full phase space for I ^ IB, 
and 

R m = 214.5MeV/c < p t < 231.5MeV/c . (1.20) 

for the Bremsstrahlung contribution. Here pi stands for the modulus of the lepton three 
momentum in the kaon rest system []. We consider constant form factors V, A and write for 
the rates and for the corresponding branching ratios 



r, 

BR/ 



A SD {M K (V±A)} Nl Xj 
Tj/T tot = N{M K (V±A)} Nl Xj 



X21) 



1 This cut has been used in || for K^-y, because this kinematical region is free from background. 
We apply it here for illustration also to the electron mode K e 2 y . 
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Table 1.1: The quantities Xi,Nj. SD ± and INT 1 * 1 are evaluated with full phase space, IB 
with restricted kinematics (|1.20|) . 





SD+ 


SD~ 


INT+ 


INT - 


IB 






1.67 • 10" 2 


1.67- lO- 2 


-8.22 • 1(T 8 


3.67- 1(T 6 


3.58- 1(T 6 






1.18 ■ 10" 2 


1.18 ■ 1(T 2 


-1.78- 10~ 3 


1.23 ■ 1(T 2 


3.68- lO" 2 






2 


2 


1 


1 








with 

N = A SD /T tot = 8.348 ■ 10~ 2 . (1.22) 

The values for Nj and Xi are listed in table |1.1| . 

To estimate Tj and BR/, we note that the form factors V, A are of order 

M K (V + A) ~ -10 -1 , M K (V - A) ~ -4 • 10~ 2 . (1.23) 

From this and from the entries in the table one concludes that for the above regions Ri, 
the interference terms INT 1 * 1 are negligible in K e2l) whereas they are important in K^j ■ 
Furthermore, IB is negligible for K e2l , because it is helicity suppressed as can be seen from 
the factor mj in A m . This term dominates however in K^y 

1.3 Determination of A(W 2 ) and V(W 2 ) 

The decay rate contains two real functions 

F ± {W 2 ) = V(W 2 ) ± A(W 2 ) (1.24) 

as the only unknowns. In Figs. ( |1.1| , |1.2|) we display contour plots for the density distributions 
fm, ■ ■ ■ , /int± f° r I — Z 1 , e - These five terms have obviously very different Dalitz plots. 
Therefore, in principle, one can determine the strength of each term by choosing a suitable 
kinematical region of observation. To pin down F^, it would be sufficient to measure at each 
photon energy the interference term INT ± . This has not yet been achieved so far, either 
because the contribution of INT ± is too small (in K e 2 7 ) , or because too few events have 
been collected (in K^). On the other hand, from a measurement of SD ± alone one can 
determine A, V only up to a fourfold ambiguity: 

SD± - {(V, A); -(V, A); (A, V); -(A, V)} . (1.25) 

In terms of the ratio 

IK = A/V (1.26) 

this ambiguity amounts to 

SD ± - { lK ; 1/ 7JC } . (1.27) 

Therefore, in order to pin down the amplitudes A and V uniquely, one must measure the 
interference terms INT 1 * 1 as well. 
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Figure 1.1: Contour plots for f lB ,. . . , / INT ± [-^27]- The numbering on the lines points towards 
increasing modulus. The normalization is arbitrary. 
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Figure 1.2: Contour plots for f lB ,. . . , / INT ± [i^ e 2 7 ]- The numbering on the lines points towards 
increasing modulus. The normalization is arbitrary. 
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1.4 Previous experiments 



K + — > e + ^ e 7 

The PDG uses data from two experiments 0, H, both of which have been sensitive 
mainly to the SD+ term in (|TT5|) . In 0, 56 events with £ 7 > 100 MeV, E e + > 236 MeV 
and 9 e + 1 > 120° have been identified, whereas the later experiment |§ has collected 51 
events with i? 7 > 48 MeV, E e + > 235 MeV and # e + 7 > 140°. In these kinematical regions, 
background from K + — > e + z/ e 7r° is absent because i?™ ax (i^ e 3) = 228 MeV. The combined 
result of both experiments is Q @ 

r(SD + )/r(ir M2 ) = (2.4 ± 0.36) • 10" 5 . (1.28) 

For SD~, the bound 

r(SD-)/r total < 1.6 • 10" 4 (1.29) 

has been obtained from a sample of electrons with energies 220 MeV < E e < 230 MeV ||. 
Using ( |Piyp2D , the result flOg) leads to 



M K | V + A |= 0.105 ± 0.008 . (1.30) 

The bound (1.29) on the other hand implies |2j 

\V-A\/\V + A\< y/Tl, (1.31) 

from where one concludes that 7^ is outside the range —1.86 to —0.54, 

7^ ^ [-1.86,-0.54] . (1.32) 

As we already mentioned, the interference terms INT 1 * 1 in K — > ez/ e 7 are small and can hardly 
ever be measured. As a result of this, the amplitudes A, V and the ratio 7^ determined from 
K e2l are subject to the ambiguities (|1.25|) , (|1.27| ). 

K + -» /i + ^7 

Here, the interference terms INT 1 * 1 are nonnegligible in appropriate regions of phase space 
(see Figs. (|l.l|J1.2j )). Therefore, this decay allows one in principle to pin down V and A. The 
PDG uses data from two experiments @, §J- In 0, the momentum spectrum of the muon 
was measured in the region (p..20[) . In total 2 ± 3.44 SD + events have been found with 216 



MeV/c < < 230 MeV/c and £ 7 > 100 MeV, which leads to 

M K I V + A |< 0.16 . (1.33) 

2 In all four experiments || ||, ](J discussed here and below , the form factors A and V have been treated 
as constants. 
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Table 1.2: Measured branching ratios T(K lvi r y)/Y t 
the K^-y data from [||, [|. The last column corresponds 



The K e2l data are from R 



to the cut (ODD 





SD + 


SD~ 


INT+ 


SD" + INT" 


total 




(1.52 ±0.23) • 10~ 5 


< 1.6 • 10~ 4 










< 3 • HT" 




< 2.7- 10~ b 
(modulus) 


< 2.6 • 1(T 4 
(modulus) 


(3.02 ±0.10) ■ 10~ 3 



In order to identify the effect of the SD terms, the region 120 MeV/c < <150 MeV/c 
was searched. Here, the background from K^ 3 decays was very serious. The authors found 



142 candidates and conclude that 



1.77 < M K (V — A) < 0.21. 



1.34) 



The result Ql.33 ) is consistent with ( 1.3C ), and the bound ( |1 .34 ) is worse than the result 
( |1.31| ) obtained from K e2l . The branching ratios which follow 
played in table [L2], where we also show the K e2l results ||. 



]] from ( |1.33|JiT34"|) are dis- 
The entry SD"+INT~ for 



is based on additional constraints from K e2l ||. 



1.5 Theory 

The amplitudes A(W 2 ) and V^H^ 2 ) have been worked out in the framework of various ap- 
proaches, viz., current algebra, PCAC, resonance exchange, dispersion relations, .... For a 
rather detailed review together with an extensive list of references up to 1976 see @]. Here, 
we concentrate on the predictions of V, A in the framework of CHPT. 

A) Chiral expansion to one loop 

The amplitudes A and V have been evaluated || || in the framework of CHPT to one 
loop. At leading order in the low-energy expansion, one has 

A = V = 0. (1.35) 

As a consequence of this, the rate is entirely given by the IB contribution at leading order. 
At the one-loop level, one finds 



8tt 2 F 

lK = 32n 2 (L r 9 + L r w ) , (1.36) 
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Table 1.3: Chiral prediction at order p 4 for the branching ratios T(K — > £^7)/T total . The cut 
used in the last column is given in Eq. ( |1.2(Jp . 





SD+ 


SD" 


INT+ 


INT- 


total 




1.30 • 10" 5 


1.95 ■ 10" 6 


6.64- 10~ 10 


-1.15 ■ lO" 8 


2.34- 10" 6 


Kfi2j 


9.24- 10~ 6 


1.38 • 10- 6 


1.44 • 10~ 5 


-3.83- 10- 5 


3.08- 10~ 3 



Lg and L"[ are the renormalized low-energy couplings evaluated at the scale \x (the 



where 

combination L 



9 + L T W is scale independent). The vector form factor stems from the Wess- 
Zumino term which enters the low-energy expansion at order p 4 , see chapter [|. 

Remarks: 



(i) At this order in the low-energy expansion, the form factors A, V do not exhibit any 

W /2 -dependence. A nontrivial jy 2 -dependence only occurs at the next order in the 
energy expansion (two-loop effect, see the discussion below). Note that the available 
analyses of experimental data of K — > lv\^ decays || 0, [|, use constant form factors 
throughout. 

(ii) Once the combination Lg + L\q has been pinned down from other processes, Eq. ( |1.36| ) 

allows one to evaluate A, V unambiguously at this order in the low-energy expansion. 
Using L 9 + L 10 = 1.4 ■ 10 -3 and F = F n , one has 

M K (A + V) = -0.097 
M K {V-A) = -0.037 

j K = 0.45 . (1.37) 

The result for the combination (A + V) agrees with ( |1.30| ) within the errors, while jk 
is consistent with (|1.32|) . 



We display in table |1.3| the branching ratios BRf (|1.21|) which follow from the prediction 
1.371 ) • These predictions satisfy of course the inequalities found from experimental data (see 



table |OD. 
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Figure 1.3: The rate dP(x)/dx in ( |1.39| ), evaluated with the form factors ( |1.38| ) and N tot = 
9 • 10 9 . The solid line corresponds to M K * = 890 MeV, M Kl = 1.3 GeV. The dashed 
line is evaluated with Mx* = 890 MeV, M%_ x = oo and the dotted line corresponds to 
Mk* = = oo. The total number of events is also indicated in each case. 

B) W^-dependence of the form factors 

The chiral prediction gives constant form factors at order p A . Terms of order p 6 have not 
yet been calculated. They would, however, generate a nontrivial W 2 - dependence both in 
V and A. In order to estimate the magnitude of these corrections, we consider one class of 
p 6 - contributions: terms which are generated by vector and axial vector resonance exchange 
with strangeness |7|, [11]], 



V(W 2 ) = Xr, 9 , A (W 2 ) = 4 9 (1-38) 

v ' l-W 2 /M K * 2 ' v ; l-W 2 /M K 2 y ' 



where V, A are given in (|1.36| ) . We now examine the effect of the denominators in ( |1.38| ) in 
the region y > 0.95, x > 0.2 which has been explored in K + — > e + v e r y @. We put m e = 
and evaluate the rate 

dP{x) N tot 



p SD +(x,y)dy (1.39) 

dx 1 tot 7^=0.95 

where N tot denotes the total number of K + decays considered, and T~l = 1.24 • 10~ 8 sec. 
The function ^j^- is displayed in Fig. ( |1.3| ) for three different values of M K * and M Kl , 
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with N tot = 9 • 10 9 . The total number of events 

N P = f 1 dP{x) (1.40) 

Jx=0.2 

is also indicated in each case. The difference between the dashed and the dotted line shows 
that the nearby singularity in the anomaly form factor influences the decay rate substantially 
at low photon energies. The effect disappears at x — *■ 1, where W 2 = M|-(l — x) —>■ 0. To 
minimize the effect of resonance exchange, the large x-region should thus be considered. The 
low x-region, on the other hand, may be used to explore the W /2 -dependence of V and of A. 
For a rather exhaustive discussion of the relevance of this W 2 - dependence for the analysis 
of Ki 2l decays we refer the reader to Ref. 0. 

1.6 Improvements at DAFNE 

Previous experiments have used various cuts in phase space in order (i) to identify the 
individual contributions IB, SD 1 * 1 , INT ± as far as possible, and (ii) to reduce the background 
from K[ 3 decays. This background has in fact forced so severe cuts that only the upper end 
of the lepton spectrum remained. 

The experimental possibilities to reduce background from Ki 3 decays are presumably 
more favourable with today's techniques. Furthermore, the annual yield of 9-10 9 K + decays at 
DAFNE is more than two orders of magnitude higher than the samples which were available 
in 0, |3|, |5], ||. This allows for a big improvement in the determination of the amplitudes A and 
V, in particular in K^-y decays. It would be very interesting to pin down the combination 
L9 + L10 of the low-energy constants which occur in the chiral representation of the amplitude 
A and to investigate the W /2 -dependence of the form factors. 
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2 The decays K± -> Pvl'+l'- 

Here we consider decays where the photon turns into a lepton-anti-lepton pair, 



K + - 


-> e 4 


"z//i H 


> 


(2.1) 


K + - 


-> ^ 


w 




(2.2) 


K + - 




4 

i/e 


e~ 


(2.3) 


K + - 


-> 




V- 


(2.4) 



2.1 Matrix elements 

We start with the processes Q2.1] ) and 



The matrix element is 
where 



tf+(p) ^ l + {pi)v{Vu)l' + {Pl)l'-{V2) 

(1,1') = (e,fi) or (/x,e). (2.5) 
T = -iGjreOp {F K r - (2.6) 



m,«(p„)(l + 75) <! — "J - ^— — -5 



with 



2//' - (f 2pf • « • ; ' 
2pq — g 2 2p ; g + q 2 

H pfl = W^qapp- A^qWgW -W p q p ) 

-A 2 (q 2 g p » - q p q p ) - A 4 (qWq p - q 2 W p )W p (2.7) 

A 2Fk F * {q2) ~ 1 I A (2Z) 

M ~ Mjt-W* ? + As - (2 ' 8) 

The form factors Ai(q 2 , W 2 ), Vi(q 2 , W 2 ) are the ones defined in appendix y. Fy(q 2 ) is the 
electromagnetic form factor of the K + . Finally the quantity e M stands for 

e" = ^u(p 2 )Yv( Pl ) , (2.9) 

and the four-momenta are 

q = Pi +P2, W = pi +p u = p - q (2.10) 

such that g M e M = 0. 

In order to obtain the matrix element for fl2.3|) and (|2.4| ), 



K + {p) - l + {pi)v{p,)l + {pi)r (pa) , (2.11) 
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one identifies mi and m\ in ( |2.6| ) and subtracts the contribution obtained from interchanging 
Pi <-> Pi ■ 

(PuPi) -> (Pi,Pl) 

W -> p-q=p v +pi . (2.12) 



2.2 Decay distributions 

The decay width is given by 



dT = 2M k (2-k) 8 ^ \ T2 \ d LiPs{p;Pi,Pu,Pi,P2) (2.13) 



spins 



and the total rate is the integral over this for the case / ^ V . For the case I = V the integral 
has to be divided by the factor 2 for two identical particles in the final state. 

We first consider the case where I ^ V and introduce the dimensionless variables 

2pq 



x 



Ml 
2piP 
Ml 

_l_ 

Ml 
mf 
Ml 



ri = ™l . (2.14) 

1 M 2 K y 1 



Then one obtains, after integrating over pi and p2 at fixed q |12 , 



dY K+ ^ l+vV+v - = a 2 G 2 F \V us \ 2 M 5 K F(z,r' l ){-J2TlT' M \dxdydz 



spins 



F(z y t) = _!_/ 1 + **U_*i 

v ' " 192tt 3 2 I z I V ^ 



T = M?{F K ir -W v l v ). (2.15) 

The quantity j— S sp ms is displayed in appendix 0. This result allows one to evaluate, 

e.g., the distribution dT/dz of produced pairs rather easily. The kinematically allowed 
region is 

Ar[ < z < 1 + n - 2^/¥i 
2\fz < x < 1 + z — n 
A-B< y <A + B (2.16) 
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with 



{2 - x){l + z + ri - x) 
2(1 + z-x) 

B = (l + *-*-WflEg . (2.17) 
2(l + z-x) V ; 

The case I = I' is slightly more elaborate. We feel that it does not make sense to display 
the term Jospins \T\ 2 because it is of considerable complexity in the general case when all 
the form factors A^, V\ and Fy are q 2 and W 2 dependent. The expression together with 
the Monte Carlo program to do the phase space integrals is available on request from the 
authors. 



2.3 Theory 

The form factors A4, V\ and Fy have been discussed in all kinds of models, Vector Meson 
Dominance, hard meson, etc.. For a discussion see Ref. f7j. We will restrict ourselves to the 
predictions in the framework of CHPT. 

To leading order we have 

Vi = 

A 1 = A 2 = A 3 = . (2.18) 

We also have Fy = 1. The rate here is entirely given by the inner Bremsstrahlung contri- 
bution. At the one-loop level several form factors get non-zero values Pfl 



1 



8tt 2 F 
4 



M = --(L r 9 + U 10 ) 

'I I ■ \ I , 

A 2 = - 
A 3 = 



2F K (F«(q 2 ) - 1) 



q 2 



F«(q 2 ) = l + H 7Tn (q 2 ) + 2H KK (q 2 ) . (2.19) 

These results obey the current algebra relation of Ref. The function Fy(q 2 ) does, 

however, deviate somewhat from the linear parametrization often used. The function H{t) 
is defined in appendix |B|. 

The fact that the form factors at next-to-leading order could be written in terms of the 
kaon electromagnetic form factor in a simple way is not true anymore at the p 6 level. The 
Lagrangian at order p 6 contains a term of the form 

tv{D a F^D^F R ^u} (2.20) 
that contributes to A 2 and A 3 but not to the kaon electromagnetic form factor, Fy(q 2 ). 
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Table 2.1: Theoretical values for the branching ratios for the decay K + — > jj, + ue + e for 
various cuts. 





tree level 


form factors as given by CHPT 


full phase space 


2.49 • 10" 5 


2.49 • 10~ 5 


z < 10" 3 


2.07- 10~ 5 


2.07- 10~ 5 


z > 10~ 3 


4.12 ■ 10~ 6 


4.20 ■ lO" 6 


z > (20 MeV/M K f 


3.15 • lO" 6 


3.23 • 10~ 6 


z > (140 MeV/M K f 


4.98 ■ 10~ 8 


8.51 • 10~ 8 


x > 40 MeV/M K 


1.58 ■ 10~ 5 


1.58 ■ 10~ 5 



2.4 Numerical results 

Using the formulas of the previous subsections and appendix [D] we have calculated the rates 
for a few cuts, including those given in the literature. For the case of unequal leptons, the 



results are given in table |2.1| for the decay K + — > fi + ue + e . These include the cuts used in 
Refs. |Lj§ and 0, x > 40 MeV/M K and z > (140 MeV/M K ) 2 , respectively. It can be seen 



that for this decay most of the branching ratio is generated at very low electron-positron 
invariant masses. As can be seen from the result for the cuts used in Ref. |Tj|, the effect of 



the structure dependent terms is most visible at high invariant electron-positron invariant 
mass. Our calculation, including the effect of the form factors agrees well with their data. 
We disagree, however, with the numerical result obtained by Ref. |L2| by about an order of 
magnitude. 

For the decay K + — > e + b>fi + {i~, we obtain for the tree level or IB contribution a branching 
ratio 

BR IB {K + -> e + vfi + fi-) = 3.06 • 10~ 12 (2.21) 
and, including the form factors, 

BR total {K + -> e+zV/O = 1.12 ■ 10" 8 . (2.22) 

Here the structure dependent terms are the leading contribution since the inner Bremsstrahlung 
contribution is helicity suppressed as can be seen from the factor mi in L M . 

For the decays with identical leptons we obtain for the muon case a branching ratio of 

BR total {K + -> /i+zV/O = 1-35 ■ 10~ 8 (2.23) 

for the full phase space including the effects of the form factors. The inner Bremsstrahlung 
or the tree level branching ratio for this decay is 

BR IB (K + -»■ n+vfi+iT) = 3.79 ■ 10" 9 . (2.24) 

For the decay with two positrons and one electron the integration over full phase space for 
the tree level results is very sensitive to the behaviour for small pair masses. We have given 
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Table 2.2: Theoretical values for the branching ratios for the decay K + — > e + z/e + e for 
various cuts. 





tree level 


form factors as given by CHPT 


full phase space 


« 4 • 10" 9 


1.8 • 10" 7 


z, z x > 10- 3 


3.0 • 10" 10 


1.22 ■ 10~ 7 


z, z x > (50 MeV/M K f 


5.2 ■ 10- 11 


8.88- 10- 8 


z, Zl > (140 MeV/M K f 


2.1 • 10" 12 


3.39- 10~ 8 



the tree level and the full prediction, including form factor effects in table [2.2| . The cuts are 
always on both invariant masses : 



z 



( Pl +p 2 ) 2 /M 

(P,+ P2 ) 2 /M: 



K 
K 



(2.25) 



The values for the masses used are those of K + and tt + . For L$ and L\q we used the values 
given in chapter |l], 



Ll(M p ) 



6.9 • 10 



-3 



-5.5 • 10" 



(2.26) 



2.5 Present experimental status 



Only decays with an electron positron pair in the final state, decays (2J!) and (|2.3j ), have 
been observed. 

Both have been measured in the same experiment |T3 . The decay K + — > yU + z/e + e" was 



measured with a branching ratio of (1.23 ± 0.32) • 10" 7 with a lower cut on the electron 
positron invariant mass of 140 MeV. The measurement is compatible with our calculation 
including the form factor effects for the relevant region of phase space. This measurement 
was then extrapolated |U| using the result of [12| to the full phase space. Since we disagree 



with that calculation, we also disagree with the extrapolation. 

In the same experiment, 4 events of the type K + — > e + ve + e~ were observed where 
both electron positron pair invariant masses were above 140 MeV . This corresponds to a 
branching ratio for this region of phase space of (2.8i 2 j) • 10~ 8 . This result is compatible 
within errors with our calculation, see table |2.2| . The matrix element of Ref. fTTH was again 



used for the extrapolation to full phase space [|13j. Apart from our numerical disagreement, 
the calculation of Ref. (TTD was for the case of non-identical leptons and cannot be applied 
here. 



For the decay K + — > ^i + v^i + ^i an upper limit of 4.1 • 10 7 exists [^4 . This upper limit 
is compatible with our theoretical result, Eq. ( |2.23| ). 
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The decay K + — > e + vji + \i~ has not been looked for so far and should be within the 
capabilities of DAFNE given the branching ratio predicted in the previous subsection. This 
decay proceeds almost entirely through the structure dependent terms and is as such a good 
test of our calculation. 



2.6 Improvements at DAFNE 

The decays discussed in this section, K + — > are complementary to the decays 



K + — > l + u^. As was the case for the analogous decay, 7r + — > e + ue + e~ [15], it may be 
possible to explore phase space more easily with this process than with K + — ► l + vy to 
resolve ambiguities in the form factors. 



As can be seen from our predictions, tables ETTI and all the decays considered in this 



section should be observable at DAFNE. Large improvements in statistics are possible since 
less severe cuts than those used in the past experiments should be possible. In the decays 
with a pair and the decay K + — > e + ve + e~ the effects of the form factors are already 

large in the total rates and should be easily visible at DAFNE. In the decay K + — ► fi + ue + e~ 
most of the total rate is for small invariant mass of the pair and is given by the inner 
Bremsstrahlung contribution. There are, however, regions of phase space where the form 
factor effects are large and DAFNE should have enough statistics to be able to study these 
regions. 
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3 K13 decays 



The decay channels considered in this section are 

K + (p) -> Ap')1 + (pM(p») (3-1) 

K°(p) - v-{p')l + (Pi>i{P») [KU (3-2) 

and their charge conjugate modes. The symbol I stands for fi or e. We do not consider 
electromagnetic corrections and correspondingly set a = throughout this section. 

3.1 Matrix elements and kinematics 

The matrix element for has the general structure 

T = ^KWV,p) (3.3) 

with 

P = u(p u )Y(l-l 5 )v( P i) 



F, + (p',p) = < 7r (p') | V;f~ i5 (0) | K + (p) > 



= ^[{p' + p),f^\t) + {p-p'),f^\t)}. (3.4) 

To obtain the matrix element for , one replaces Fj~ by 

F,°(p',p) = <n-(p')\V^(0)\K°(p)> 

= (p' + p),ff 7T ~(t) + (p-p>),f K °--(t). (3.5) 



The processes fl3.1|) and (|3.2|) thus involve the four form factors f± +7T °(t), f±° w (t) which 
depend on 

t={p' -pf = {p l+ p y )\ (3.6) 

the square of the four momentum transfer to the leptons. 

Let f± n = f± +7T ° or f± ° n . is referred to as the vector form factor, because it 

specifies the P-wave projection of the crossed channel matrix elements < | V^~ l5 (0) \ 
K + , 7r° in >. The S'-wave projection is described by the scalar form factor 



ft(t) = f«*(t) + M ^_ Ml f-\t) ■ (3-7) 



Analyses of data frequently assume a linear dependence 

t 



/5(f) = /^(o) 



1 + A+,o 



M 2 



(3.8) 
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For a discussion of the validity of this approximation see [JT6|, [I[ and references cited therein. 
Eq. O leads to a constant f^{t) , 

M 2 — M 2 

f Kn (t) = f^(0) = /^(0)(A - A + ) K \ (3.9) 

The form factors f±$(t) are analytic functions in the complex t-plane cut along the positive 
real axis. The cut starts at t = (M K + M n ) 2 . In our phase convention, the form factors are 
real in the physical region 

mf < t < (M K - M w ) 2 . (3.10) 
The kinematics of (spin averaged) decays needs two variables, for which we choose 

y = 2p Pl /M 2 K , z = 2pp'/M 2 K = (-t + M 2 + M 2 K )/M 2 K . (3.11) 

In the K rest frame, y (z) is proportional to the charged lepton (pion) energy, 

y = 2E l /M K , z = 2EJM k . (3.12) 

The physical region for y and z is 

2^/ri < y < 1 + n - r w 
A(y)-B(y)< z <A(y) + B(y) 

A(y) = ( 2 - y)(l + n + ry - y)/[2(l + n - y)] 

B{y) = ^y 2 - 4r,(l + r, - ry - y)/[2(l + n - y)} 

n = m 2 /M 2 K ,r n = M 2 /M 2 K . (3.13) 

or, equivalently, 

2^fr\ < z < 1 + ry - n 
C{z)-D{z)< y <C(z)+D(z) 

C{z) = {2-z){l + r w + n-z)/[2(l + r^-z)} 

z)] . 

(3.14) 



D{z) = Jz 2 -4r 7T (l + r 7T -r l -z)/[2(l + r 7r -z)} 



3.2 Decay rates 

The differential decay rate for is given by 



dT = 2Mk{2tt) 5 ^ \ T \ 2(1 lips{p;Pi,Pu,p') • (3.15) 



spins 

The Dalitz plot density 



OM = £r. = E W (3.16) 



spins 
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is a Lorentz invariant function which contains f± +n ° in the following form, 



MlGl I V m I 2 



P(y,z) = ' [A(f^y + Bf^°f^° + C(f^y\ (3.17) 



with 



A(y,z) = 4(z + y-l)(l-y)+r l [4y + 3z-3]-4r n + r l (r„-r l ) 
B(y,z) = 2r ; (3 - 2y - z + n - r n ) 

C{y,z) = nil + r^-z-n). (3.18) 
The quantities (A, B, C) are related to the ones quoted by the PDG |L| by 

(A,B,C) = ^(A,B,C) mG . (3.19) 

To obtain the rate for one replaces in ( |3.17| ) f± +n ° by \^2f^ 0jT . 

For convenience we also display the K^ 3 /K e3 rates evaluated in the approximation ( |3.8|) 
for the form factors, 

+ + 0.645 + 2.087A+ + 1.464A + 3.375A^ + 2.573A 2 , 

( ^' { e3) ~ 1 + 3.457A+ + 4.783A 2 , 

0.645 + 2.086A+ + 1.459A + 3.369A 2 , + 2.560A 2 
m^/T(Kj = i + 3.456A + + 4.776A^ ' (3 - 20) 

We have used the physical masses U in evaluating these ratios and M n + to scale the slope 
in both cases. The terms linear and quadratic in Ao are proportional to mf and therefore 
strongly suppressed in the electron case. We do not include them in the denominators, 
because these coefficients are smaller than 10~ 4 . The interference term AoA + is absent by 
angular momentum conservation. Furthermore, one has 

*^-) = {ao606 @ . < 3 - 21 > 

3.3 Determination of the Ki% form factors 



Measurements of the Dalitz plot distribution ( |3.17| ) of data allow one in principle to pin 



down the form factors (up to a sign) in the range m 2 < t < (M K - M n ) 2 . Measuring the 
Kps/Kes branching ratio and then using ( |3.20| ) gives a relationship between A + and Ao which 



is valid in the approximation ( |3.8|) . Furthermore, muon polarization experiments measure 
the weighted average of the ratio /^ 7r (t)//^' 7r (t) over the t range of the experiment [|I|, |T7| |. 
On the other hand, the electron modes K e3 are sensitive to only, because the other 
contributions are suppressed by the factor (m e /M^) 2 ~ 10~ 6 , see eqs. ( |3.17|) , (|3.18|) . 

Isospin breaking effects in f+ +7r °(0) and f+° n (0) play a central role in the determina- 
tion of the Kobayashi-Maskawa matrix element V us from K e3 data, see [|ll| for a detailed 



discussion of this point. In the following we concentrate on the measurement of the slopes 

^+,0- 
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3.4 Previous measurements 



We refer the reader to the 1982 version of the PDG Jl9| f] for a critical discussion of the wealth 
xperimental info] 

K e3 -exp eriment s 



of experimental information on A+ q. Here we content ourselves with a short summary. 



The A + values obtained are fairly consistent. The average values are 



*3 

K° 



K 
A + 



0.028 ±0.004 Ref.[l] 
0.030 ±0.0016 Ref.|HH 



(3.22) 



-K" M 3-experiments 

The result by Donaldson et al. |20 



A+ = 0.030 ±0.003 
A = 0.019 ±0.004 



(3.23) 



dominates the statistics in the K® 3 case. The A + value ( |3.23| ) is consistent with the K e3 



value (|3.22|) . However, the situation concerning the slope Ao is rather unsatisfactory, as the 
following (chronological) list illustrates. 



0.0341 

0.050 

0.039 

0.047 

0.025 

0.019 



± 
± 
± 
± 
± 
± 



0.0067 

0.008 

0.010 

0.009 

0.019 

0.004 



[21] 

in 

H 
m 



(3.24) 



The x 2 fit to the if° 3 data yields A+ = 0.034 ± 0.005, A = 0.025 ± 0.006 with a x 2 /DF 
88/16 [O, p. 76]! The situation in the charged mode K% is slightly better |Tj5 . 



3.5 Theory 



The theoretical prediction of Ki 3 form factors has a long history, starting in the sixties with 
the current algebra evaluation of f± +7T . For an early review of the subject and for references 



to work prior to CHPT evaluations of f± we refer the reader to p6| (see also Ref.[[27]]). Here 



we concentrate on the evaluation of the form factors in the framework of CHPT. We restrict 
our consideration to the isospin symmetry limit m u = m^, as a result of which one has 



J±,0 



(t) 



rK+; 
J±,0 



(t) = f±,o{ t ) ; m u = m d 



(3.25) 



A) Chiral prediction at one-loop order 

3 Please note that the most recent measurements of A+.o go back to 1981 [0! 
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In Ref. [fTE| , the vector current matrix elements < M' \ q^—q \ M > have been 
calculated up to and including terms of order t = {p' — p) 2 and of order m u ,md and m s in 
the invariant form factors. For reasons which will become evident below, we consider here, 
in addition to the Ki 3 form factors, also the electromagnetic form factor of the pion 

< ^ + {P') I Km(0) I * + (p) >= (p'+pYF*(t). (3.26) 
The low-energy representation for Fy(t) [|1(| ^| and /+(£) [HJ reads 



F$(t) = l + 2H mr (t) +H KK (t) 

f + (t) = 1 + ^H Kn (t) + ^H Kv (t). (3.27) 

The quantity H(t) is a loop function displayed in appendix [TJ. It contains the low-energy 
constant L$. The indices attached to H(t) denote the masses running in the loop. 

Since L 9 is the only unknown occurring in Fy(t) and in f+(t), we need experimental 
information on the slope of one of these two form factors to obtain a parameter-free low- 
energy representation of the other. 

The analogous low-energy representation of the scalar form factor is 
/o(*) = l + ^Ut-2V Kv -3^Jj Kir (t) 

+ T^f^-lV (3-28) 



A-Kn V F n 

The function J(t) is listed in appendix [B], and and A Kn stand for 



E Kw = M 2 K + Ml 

A Kn = M 2 K -Ml . (3.29) 



The measured value |18| Fk/F w = 1.22 ± 0.01 may be used to obtain a parameter-free 



prediction of the scalar form factor fo(t). 

B) Momentum dependence of the vector form factor 



In the spacelike interval < 350 MeV the low-energy representation (|3.27|) for the 
electromagnetic form factor Fy(t) is very well approximated by the first two terms in the 
Taylor series expansion around t — 0, 

F v (t) = l + ^<r 2 > v t + --- . (3.30) 
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Figure 3.1: The vector and scalar form factors /+(£) and fo(t). 



Likewise, the linear approximation 



f + (t) = f + (0)\l + -<r 2 >^t + 



(3.31; 



reproduces the low-energy representation ( 3.27 ) very well, see Fig. This is in agreement 
with the observed Dalitz plot distribution, which is consistent with a form factor linear in t. 
The charge radii are 



2 IT 12Lg 

< r 2 >y~ 



< r > 



2 ^ Kir 



F 2 32vr 2 F 2 

1 



M 2 
2 In — f + In 



M 



K 



+ 3 



2 ^ 7T 



< r > 



v 



64n 2 F 2 



3h 



Ml, 



+ 3/ii 



Ml, 



+ - In 



M 



K 



+ -In 



M 2 



where 



2 ~ Ml ' 2 " M 2 K 



, 1 (x 3 - 3x 2 - 3s + 1) , l/.x + l^ 2 

Mx) = 2 ^Tp ln " + 2 l^T, 



1 

3 



(3.32) 



(3.33) 



To evaluate these relations numerically, we use the measured charge radius of the pion: 

< r 2 >l= 0.439 ± 0.008fm 2 Ref.[29] (3.34) 
as input and obtain the prediction 

A + = -M 2 + < r 2 >* 7T = 0.031 (3.35) 
6 
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in agreement with the experimental results ( |3.22| ), ( |3.23| ) Q . From this (and from the 



considerably more detailed discussion in Ref. jL6|), one concludes, in agreement with other 
theoretical investigations [ 3"0|] , that the measured charge radii < r 2 >y and < r 2 >y n are 



consistent with the low-energy prediction. 

C) Momentum dependence of fo(t). Dashen-Weinstein and Callan-Treiman 
relations 



In the physical region of decay the low-energy representation ( |3.28 ) for the scalar 
form factor is approximated by the linear formula 

hit) = f + (0) |l + ~ < t 2 >*«t + ■ ■ j (3.36) 



to within an accuracy of 1 %. (See Fig. |3.1| ). The curvature generated by higher order 
terms is also expected to be negligible in the physical region of the decay JM. For the slope 
< r 2 >g n one obtains 



62 ~ -T92^ 2 \ 15h2 {M^) + M 2 K + M 2 k2 \Ml, 



(3.37) 



where 



, , s 3 fl + x\ 2 3x(l + x) 

h2\x) = - + —, —ma;, 

K ' 2\l-xJ (1 -x) 3 

h 2 (x) = h 2 (-} , h 2 (l) = 1, 

rh = {m u + m d )/2 . (3.38) 



This (parameter-free) prediction is a modified version of the Dashen-Weinstein relation |JT 
which results if the nonanalytic contribution 5 2 is dropped. Dashen, Li, Pagels and Weinstein 



were the first to point out that the low-energy singularities generated by the Goldstone 
bosons affect this relation. The modified relation is formulated as a prediction for the slope 
of fo(t) at the unphysical point t\ = M\ + M 2 . Their expression for this slope however has 
two shortcomings: (i) it does not account for all corrections of order M.: (ii) The slope at t\ 
differs substantially from the slope in the physical region of the decay 0, |33|, see Fig. |~2|. 

Algebraically, the correction 6 2 is of the same order in the low-energy expansion as the 
term involving F^/F^ — 1. Numerically, the correction is however small: 5 2 reduces the 

4 We do not quote an error for the result (|3.35| ), because one should estimate higher order chiral corrections 
for this purpose. 
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Figure 3.2: The normalized slopes of the vector and the scalar form factors. Curve 1: the 
normalized slope M%+df + (t)/dt. Curve 2: the normalized slope M%+dfo(t)/dt. Near the nK 
threshold t = (M K + M w ) 2 , the vector form factor behaves as /+(£) = f+(t ) + 0[(t — to)], 
whereas f (t) = fo(t ) +0[(y/t — to)]. The slope of the scalar form factor is therefore singular 
at t = (M K + Mir) 2 . 



prediction by 11 %. With Fk/F^ = 1.22 ± 0.01 the low-energy theorem (|3.37| ) implies 



< r" > 



Kit 



An 



0.20 ± 0.05fm 2 

\mI+ < r 2 >f w = 0.017 ± 0.004 
6 



(3.39) 



where the error is an estimate of the uncertainties due to higher order contributions. The 
prediction ( |3 . 3 9| ) is in agreement with the high-statistics experiment |2^] quoted in ( |3.23p 
but in flat disagreement with some of the more recent data listed in ( ft.24Q . 

In the formulation of Dashen and Weinstein [|3ll , the Callan-Treiman relation |04l states 
that the scalar form factor evaluated at t 
order m u ,m d : the quantity 

Fk 
F w 



Mk — M 2 differs from F^ / F^ only by terms of 



A CT = f (M 2 K - Mi 
is of order rh. Indeed, the low-energy representation (p.28 ) leads to 

Ml 



(3.40) 



p f _ 1 ^1 

2 ^ [J K ,{M 2 K - Ml) + -J K ,(M 2 K - M 2 )j + 0(mm s ) . (3.41) 

—3.5 • 10~ 3 . The Callan-Treiman relation should therefore hold to a 
very high degree of accuracy. If the form factor is linear from t = to t = M\ — M 2 then 



Numerically, A c 
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Table 3.1: Rates of Ki^ decays. The number of events in the third column corresponds to 
those data which are of relevance for the determination of the slope Ao of the scalar form 
factor. 





t) events 








branching 


Particle Data 


DAFNE 


improve- 






ratio 


Group 


1 year 


ment 


K + - 




3.18 • 10- 2 


10 5 


3- 10 s 


3- 10 3 


K L - 




27- 1(T 2 


4- 10 6 


3- 10 s 


70 



the slope must be very close to 



\CT 
A 



M| - Ml 



Fk 



0.019, 



(3.42) 



in agreement with ( 3.39|) and with the experimental result of Ref 
with, e.g., the value found in Ref. J2"2"| . We see no way to reconcile the value Ao 
with chiral symmetry. 



20] , but in disagreement 
0.050 



3.6 Improvements at DAFNE 

DAFNE provides the opportunity to improve our knowledge of Ki 3 decays in a very substan- 
tial manner - in particular, it should be possible to clarify the issue of the slope Ao of the 
scalar form factor /q. To illustrate, we compare in table |3.1| the hitherto obtained number of 
events (third column) with the expected ones at DAFNE (fourth column). The last column 
displays the remarkable increase in statistics obtainable at DAFNE. 
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4 Radiative decays 



The decay channels considered in this section are 

and the charge conjugate modes. We only consider real photons (g 2 = 0). 
4.1 Matrix elements 

The matrix element for has the general structure 

T = ^eV: s e»(qyl [ (V+-A;MP»)Y(l-lMPi) (4-1) 

+ I^u(p u )Y(i - 75 )K- i> x - 4)1 Avi)] = ^K- 

2p t q J " 

The diagram of Fig. ^4].a corresponding to the first part of Eq. ( |4. 1|) includes Bremsstrahlung 
off the K + . The lepton Bremsstrahlung diagram of Fig. |4.1| .b is represented by the second 
part of Eq. ( ^4.1| ). The hadronic tensors V^, A+ v are defined as 

C = * J d 4 xe^(7r°(p') | T{V^(x)lt~ i5 (0)} | K+(p)), I = V,A. (4.2) 
is the matrix element 

K = (Ap') I K 4 ^ 5 (0) I K + {p)). (4.3) 
The tensors V^ v and A* satisfy the Ward identities 

<7%t = F+ (4.4) 
<fA% = 

leading in turn to 

?MJ = , (4.5) 

as is required by gauge invariance. 

For Kf 3 , one obtains the corresponding amplitudes and hadronic tensors by making the 
replacements 



- A% (4.6) 
4° 



K+ - 




7T°- 


V + - 

fJ,V 




A + 


F + - 

V 




A + 
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To make the infrared behaviour transparent, it is convenient to separate the tensors 
V° v into two parts: 

V+ = V+ + ~ F » ( 4 - 7 ) 

p h> pq 

yO = V° + 

Due to Low's theorem, the amplitudes V^/° are finite for q — > 0. The axial amplitudes 

are automatically infrared finite. The Ward identity ( |4.4|) implies that the vector amplitudes 

V^7° are transverse: 

?V+f> = 0. (4.8) 

For on-shell photons, Lorentz and parity invariance together with gauge invariance al- 
low the general decomposition (dropping the superscripts +,0 and terms that vanish upon 
contraction with the photon polarization vector) 

V = «W(VY + A 2 qW)+te, Xpa p'yW°(A 3 W u + A 4 p' u ) 
F v = C lP ' u + C 2 (p-p') u 

W = Pl+Pu- 



With the decomposition (|4.7| ) we can write the matrix element for in ( |4.1|) in a form 
analogous to Eq. (|1.2|) for : 

T = ^eV: s e^qyl(VX-A + Ju(p u )r(l-l 5 HPi) (4-10) 



pq 2piq 



v{pi 



The four invariant vector amplitudes V\, . . . , V4 and the four axial amplitudes A\, . . . , A4 
are functions of three scalar variables. A convenient choice for these variables is 



£ 7 = pq/M K , E n = pp'/M K , W = Vw* (4.11) 

where W is the invariant mass of the lepton pair. The amplitudes C\, C 2 can be expressed in 
terms of the K\ 3 form factors and depend only on the variable (p—p') 2 = M] c +M 2 —2MkE- !T . 
For the full kinematics of Ki 3y two more variables are needed, e.g. 

Ei = ppi/M K , x = piq/M 2 K . (4.12) 
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Figure 4.1: Diagrammatic representation of the ifjjL amplitude. 

The variable x is related to the angle 8^ between the photon and the charged lepton in the 
K rest frame: 

xM 2 K = E y (Ei - y/Ef-mf cos 9 h ). (4.13) 

T invariance implies that the vector amplitudes V\, . . . , V4, the axial amplitudes A\, . . . , A4 
and the K13 form factors C\,C2 are (separately) relatively real in the physical region. We 
choose the standard phase convention in which all amplitudes are real. 

For 6^ — > (collinear lepton and photon), there is a lepton mass singularity in ( [4.1Q 
which is numerically relevant for I = e. The region of small E 7 , 9^ is dominated by the 
K13 matrix elements. The new theoretical information of Ki^ decays resides in the tensor 
amplitudes V^ v and A^ v . The relative importance of these contributions can be enhanced 
by cutting away the region of low E 7 , 9i T It may turn out to be of advantage to reduce the 
statistics by applying more severe cuts than necessary from a purely experimental point of 
view. 



4.2 Decay rates 

The total decay rate is given by 

T(K - vrW) = 2M * / d LIPS (p;p',p h p„,q) £ | T | 2 (4.14) 



spins 



in terms of the amplitude T in ( |4.1|) . The square of the matrix element, summed over photon 
and lepton polarizations, is a bilinear form in the invariant amplitudes Vi, . . . , V4, Ai, . . . , A4, 
Ci, C%. Pulling out common factors, we write ( 4.14Q in the form 

4-n>C 2 IV I 2 r 

T(K - vr/z/7) = ^y^J J d LI p S (p;p',Pi,P„q) SM. (4.15) 
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SM, the reduced matrix element squared, is given in App. [E] as a function of scalar products 
and invariant amplitudes. For the actual numerical calculations, we have found it useful to 
employ a tensor decomposition different from the one in Eqs. ( [4. 7| ) and ( |4.9| ) 



V, 



+B 6 W»W V + B^ W V + Brp'^ . 



(4.16) 



One advantage is that ( |4.16| ) applies equally well to both charge modes while (|4.7| ) does not. 
Moreover, the expression for SM in App. [E] is more compact when written in terms of the 
Bi. In the numerical evaluation of the amplitudes, gauge invariance can of course be used 
to express three of the Bi in terms of the remaining ones and of C\,C%. 

To get some feeling for the magnitude of the various decay rates, let us first consider 
the tree level amplitudes to lowest order p 2 in CHPT. With the sign conventions of Ref.|35l 
exhibited in chapter [l], these amplitudes are || ^] : 



K + 



K° 



V. 



A + 



1 

V2 


1 

72 



(p' + W)^2p' + W) u 



pq 



(4.17) 



(P + P')u 



F?. 



+ 



(P + P') 



p'^{2p' + 2q + W\ 
p'q 



(4.18) 



In table |4. 1| the correspondin 
E y > 30MeV and 9 h > 20 



branching ratios are presented for the four decay modes for 

7T ± Z =F Z/7). 



For K® 3 , the rates are to be understood as T(Kl 



The number of events correspond to the design values for DAFNE (cf. App. [A] 



4.3 Previous experiments 

The data sample for i*Q 37 decays is very limited and it is obvious that DAFNE will be able to 



make significant improvements. The present experimental status is summarized in table [O 



A comparison between tables ^4.1| and |4.2| shows the tremendous improvement in statistics 
to be expected at DAFNE. We shall come back to the question whether this improvement 
will be sufficient to test the standard model at the next-to-leading order, 0(p 4 ), in CHPT. 
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Table 4.1: Branching ratios for tree level amplitudes for E 1 > 30MeV and 6^ > 20° in the 
K rest frame. 



decay 


BR(tree) 


#events/yr 




2.8 x 10" 4 

1.9 x 10~ 5 
3.6 x 1(T 3 
5.2 x 10" 4 


2.5 x 10 6 
1.7 x 10 5 
4.0 x 10 6 
5.7 x 10 5 



Table 4.2: Experimental results for i^ 37 decays 



decay 



-"■£37 



exp. 



E, 



7, rain 



# events 



BR 



(2.7 ±0.2) x 1(T 4 
(3.7 ± 1.3) x 10~ 4 
(2.3 ± 1.0) x 10~ 4 
< 6.1 x 10~ 5 
(1.3 ±0.8) x IP' 2 



m 

38] 
39 1 



39 

HI 



lOMeV 
WMeV 
SOMe^ 
30MeV 
IhMeV 



192 
13 
16 


10 



0.6 < cosfl e7 < 0.9 



COS# e7 < 0.9 



90% c.l. 



4.4 Theory 

Prior to CHPT, the most detailed calculations of K^-y amplitudes were performed by Fearing, 
Fischbach and Smith [^TJ using current algebra techniques. 

In the framework of CHPT, the amplitudes are given by ( [4.1 7| ) and ( [4.18 ) to leading 
order in the chiral expansion. 



A) CHPT to 0(p 4 



There are in general three types of contributions [ 55] : anomaly, local contributions due 
to £ 4 and loop amplitudes. 

The anomaly contributes to the axial amplitudes 



.4 



[IV 



A 



[XV 



/y2 f w<f(V±wr± 1 



16tt 2 F 2 

i 

_ 8tt 2 F 2 



e u x ar W v p ,x q p W 



(4.19) 



'/iv pa 



q p W a . 



The loop diagrams for Ki^ are shown in Fig. [4.2| . We first write the matrix element 
in terms of three functions ftif^i ft which will also appear in the invariant amplitudes 
. Including the contributions from the low-energy constants L 5 , L 9 in £ 4 , the K i3 matrix 
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Figure 4.2: Loop diagrams (without tadpoles) for K l3 at 0(p 4 ). For K^ p the photon must 
be appended on all charged lines and on all vertices. 



element Fj~ is given by 
F + 



A, 



4L 9 



1 



(M 2 K -M*-t)f 2 + (t) + ti(t) (p-p% 



V2F 2 t ^ 



4) 



i=i 

3 



2Bi(t) 



(t + A^AjJjjt) 
2t 



^500 - 



256tt s 
Mj -mj , t 



M n M%M v 
= (P-P') 2 ■ 



ciAjJjit) 



(4.20) 



Lg is a scale independent coupling constant and we have traded the tadpole contribution 
together with L5 for Fk/F v in fs(t). The sum over I corresponds to the three loop diagrams 
of Fig. [12] with coefficients c[, c^, C3 displayed in table |Q[ We use the Gell-Mann-Okubo 
mass formula throughout to express Ml in terms of M&,M%. The functions J/(t) and B T 2 {t) 
can be found in App. IB]. 
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Table 4.3: Coefficients for the loop amplitudes corresponding to the diagrams I = 1, 2, 3 
in Fig. |4~2] . All coefficients c\ must be divided by 6^/2F 2 . 



I 


Mr 


mi 


r 1 


r 1 


r 1 


1 




M n 


1 


-2 


-M£ - 2Ml 


2 




M v 


3 


-6 


-M\ - 2W- 


3 




M K 





-6 


-6M 2 



The standard form factors f-{t) as given in the previous section 16 are 



/+(*) = ^fiW 



(4.21] 



/-(*) - T2 m 



l-Ml 



t)f+(t) + 2f+(t)-f+(t) 



It remains to calculate the infrared finite tensor amplitude V£ v . The invariant amplitudes 
can be expressed in terms of the previously defined functions and of additional 
amplitudes Ii,I 2 ,l3. Diagrammatically, the latter amplitudes arise from those diagrams 
in Fig. ^| where the photon is not appended on the incoming K + (non-Bremsstrahlung 
diagrams). The final expressions are 



K 4 



V 3 + 



h+ P 'W q f+(W*) + f+(W* 



pq 
pq 



p'w q f+(w;) + f+(w t 



p'Wf+(W 2 ) + f+(W 2 ) -p'W q f.f(W q 2 ) - f+(W c 

- fHK) 



pq 

W + q = p — p' 



(4.22) 



The amplitudes h, h, h in Eq.( |4.22| ) are given by 
h 



AqW :(L 9 + LT ) + ^L 9 



V2F 2 



V2F 2 



+ 



E [« + A i)(cl + 4) - Mp'w q + c{) 



i=i 



(4 - 4) 



2(4 - c[) 



(W 2 - A/) J/ 



2W 2 



2Gj) 



p'W.iW^-A 2 )^ 



W 2 

q 



w 2 

q 



+ AB I 2 )+p'(W-q)L I m 



qW 



p'giFj - (W 2 + A / )G / ) + p'W{Bi - Bi)] } 
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Table 4.4: Coefficients for the loop amplitudes corresponding to the diagrams 1 = 1,2,3 
in Fig. |4~2] . All coefficients c\ must be divided by 6v^2-F 2 - 



/ 


M 7 


mj 


r 1 


r 1 




1 


M K 







-3 


-3M| 


2 


M K 


M v 


6 


-3 


M| + 2M 2 


3 


M n 


M K 


4 


-2 


2 Mj s + 2.U: 



h 
h 

Fj 
Gi 



8U 



y^F 2 qW 

3 



Fj - {W 2 + A/)G/ 



i;o(/l)+ ^ ln w. 



1 l^ 2 



(4.23) 



256t^ 



2 



S/ 1 

327r 2 A 7 n Mf 

^2 



m; 



^C(iy„ 2 , W 2 , M 2 , m 2 ) + — — {W 2 + A/) J/ - {W 2 + Aj) J, 



+ 



64vr 2 



5 2 7 = ^(W^ 2 ), B\ = Bl(W 2 ). 



The function C(W 2 ,W 2 , Mj,m 2 ) is given in App. [B|. All the invariant amplitudes 
Vi, . . . , \4 + are real in the physical region. Of course, the same is true for the matrix 
element Fj~. 

The amplitude has a very similar structure. Both the Kf z matrix element F® and 
the infrared finite vector amplitude V® can be obtained from the corresponding quantities 
F+ and V+ v by the following steps: 

• interchange p' and —p ; 

• replace — by — in f£; 



insert the appropriate coefficients c\ for Kf 3 listed in table [O 
multiply and by a factor — \/2. 



B) Numerical results 
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Table 4.5: Branching ratios and expected number of events at DAFNE for K { 





BR 


^events/yr 


full 0(p 4 ) amplitude 
tree level 

0{p A ) without loops 


3.0 x 1(T 4 
2.8 x 10" 4 
3.2 x 10" 4 


2.7 x 10 6 
2.5 x 10 6 
2.9 x 10 6 






BR 


#events/yr 


full 0(p 4 ) amplitude 
tree level 

0{p A ) without loops 


2.0 x 10" 5 
1.9 x lO" 5 

2.1 x 10" 5 


1.8 x 10 5 
1.7 x 10 5 

1.9 x 10 5 



Table 4.6: Branching ratios and expected number of events at DAFNE for K\ 





BR 


#events/yr 


full 0{p A ) amplitude 
tree level 

0(p A ) without loops 


3.8 x 10~ 3 
3.6 x 10~ 3 
4.0 x 10~ 3 


4.2 x 10 6 
4.0 x 10 6 
4.4 x 10 6 






BR 


^events/yr 


full 0(p 4 ) amplitude 
tree level 

0{p A ) without loops 


5.6 x 10" 4 
5.2 x 10~ 4 
5.9 x 10" 4 


6.1 x 10 5 
5.7 x 10 5 
6.5 x 10 5 



In calculating the rates with the complete amplitudes of the previous subsection, we use 
the same cuts as for the tree level rates in Sect. 



E 1 > 30MeV (4.24) 

e h > 20°. 

The physical values of M n and Mk are used in the amplitudes. M v is calculated from the 
Gell-Mann-Okubo mass formula. The values of the other parameters can be found in chapter 
HI and in appendix [A]. 

The results for and i£|L are displayed in tables and fO|, respectively. For 



comparison, the tree level branching ratios of table |4.1| and the rates for the amplitudes 
without the loop contributions are also shown. The separation between loop and counter- 
term contributions is of course scale dependent. This scale dependence is absorbed in the 
scale invariant constants L 9 ,L 10 defined in Eqs.( [4.2d|) , (|4.23|) . In other words, the entries in 
tables [OL WM for the amplitudes without loops correspond to setting all coefficients c{ in 



tables 4.3, PO equal to zero. 
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4.5 Improvements at DAFNE 



The numerical results given above demonstrate very clearly that the non-trivial CHPT effects 
of 0(p 4 ) can be detected at DAFNE in all four channels without any problem of statistics. 
Of course, the rates are bigger for the electronic modes. On the other hand, the relative size 
of the structure dependent terms is somewhat bigger in the muonic channels (around 8% 
for the chosen cuts). We observe that there is negative interference between the loop and 
counterterm amplitudes. 

The sensitivity to the counterterm coupling constants Lg, L\q and to the chiral anomaly 
can be expressed as the difference in the number of events between the tree level and the 
0(p A ) amplitudes (without loops). In the optimal case of K® 3l , this amounts to more than 
4 x 10 5 events/yr at DAFNE. Almost all of this difference is due to Lg. It will be very 
difficult to extract the coupling constant L 10 from the total rates. A more detailed study is 
needed to determine whether L w can be extracted from differential distributions. 

The chiral anomaly is more important for , but even there it influences the total 
rates rather little. Once again, a dedicated study of differential rates is necessary to locate 
the chiral anomaly, if possible at all. 

On the other hand, taking into account that Lg is already known to good accuracy (see 
chapter p]), Ki 3l decays will certainly allow for precise and unambiguous tests of the one- loop 
effects in CHPT g. 
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5 K14 decays 



In this section we discuss the decays 

K +(p) _> vr +( pi ) n -(p 2 ) i+(p t ) Vl {p v ) (5.1) 
K+(p) - 7T°(pi) 7r°(p 2 ) Z + (p,) (5.2) 

tf°(p) _> TT ^) 7T-( P2 ) Z+( Pi ) ^(P^) (5-3) 

and their charge conjugate modes. The letter Z stands for e or /i. We do not consider isospin 
violating contributions and correspondingly set m u = md, a. = 0. 

5.1 Kinematics 

We start with the process ( |5.1| ). The full kinematics of this decay requires five variables. We 
will use the ones introduced by Cabibbo and Maksymowicz ||42|| . It is convenient to consider 
three reference frames, namely the K + rest system (Ert), the n + ir~ center-of-mass system 
(S 27r ) and the l + ui center-of-mass system (£/„). Then the variables are 

1. s n , the effective mass squared of the dipion system, 

2. si, the effective mass squared of the dilepton system, 

3. the angle of the ir + in £2^ with respect to the dipion line of flight in £#, 

4. 81, the angle of the Z + in with respect to the dilepton line of flight in Eg", and 

5. 0, the angle between the plane formed by the pions in £# and the corresponding plane 
formed by the dileptons. 



The angles 9 n , 61 and 4> are displayed in Fig. |5.1| . In order to specify these variables more 
precisely, let pi be the three-momentum of the ir + in £2^ and pi the three-momentum of the 
Z + in £;^. Furthermore, let v be a unit vector along the direction of flight of the dipion in 
Eg-, and c(cZ) a unit vector along the projection of pi (pi) perpendicular to v(—v), 

c = (pi - w -pi)/[pi 2 - (pi • v) 2 ] 1/2 
d = (pi- vv-pi)/]p 2 - (pi ■ w) 2 ] 1/2 • 



The vectors v, c and d are indicated in Fig. 5.1. Then, one has 



Sn = (P1+P2) 2 , Si = (pi+p u f 

cos 6^ = v-pi/\pi\, cos 0i = -v ■ pi/ I pi I 

cos0 = c-d, sin = (c x tf) • d. (5.4) 
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Figure 5.1: Kinematic variables for K i4 decays. The angle 9 W is defined in E 27r , 9i in T,^ and 
in Hk- 



The range of the variables is 

4M 2 < s n <(M K - mi ) 2 

m 2 < s t < (M K - 

< 9 n , 01 < 7T, < < 2?r. 

It is useful to furthermore introduce the following combinations of four vectors 

P = Pi + P2, Q =Pi~P2 , L = pi+p„ , N = Pl - Pu 
together with the corresponding Lorentz invariant scalar products 

P 2 = s n , Q 2 = 4M 2 -s w , L 2 = s h N 2 = 2m 2 - Sl , 
PQ = 0, 

PL 



with 



1 



K °tt 



Si), 



PN 
QL 
QN 

LN 

< LNPQ > 



ziPL+ (1 - Z[ )X cos 9i, 
a v X cos6 n , 

ziQL + a w (l — z{) [PL cos9 w cosOi 
— (s^s^) 1 / 2 sin 9 n sin Q x cos <p 



m. 



X 



— (s7rS/) 1 ^ 2 o"7r(l — zi)X sin 9 n sin 9\ sinc/> 



((PL) 2 - s^) 1 ' 2 = ^ 2 (M 2 K ,s n , Sl ) 
m 2 / Sl 

(1-4M 2 / S?r ) 1 / 2 . 



(5.5) 
(5.6) 



(5.7) 



(5.8) 
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Below we will also use the variables 



t = (p 1 - pf 

u = (p 2 -p) 2 . (5.9) 



These are related to s n , S[ and 9^ by 



t + u = 2M 2 + M\ + Sl - s n 

t-u = -2ct w Xcos^ . (5.10) 



5.2 Matrix elements 

The matrix element for K + — > -k + -k^I + vi is 



where 



T = 9£v: s u(Pu)l,(l - iMPiW - ^) (5-11) 



J M = < 7r + ( Pl )^"(p 2 )out | /^ & (0) | K + (p) >; I = V,A 

= -t-^-[P^F + Q^G + L u R] (5.12) 



and eoi23 = 1- The matrix elements for the other channels ( |5.2| , |5.3| ) may be obtained from 
(|5.11| , |5TT2"|) by isospin symmetry, see below. 



The form factors F,G,R and H are real analytic functions of the three variables P1P2, 
Pip and p 2 p. Below, we will use instead the variables {s^s;,^} or {s n ,t,u}. 



Remark: In order to agree with the notation used by Pais and Treiman El and by 



Rosselet et al. jPJ, we have changed our previous convention | 45| , |46| in the definition of the 



anomaly form factor H. See also the comments after Eq. ( [5.211 ). 



5.3 Decay rates 

The partial decay rate for (|5.1|) is given by 



K \ ) spins 

The quantity J2 S pins I T | 2 is a Lorentz invariant quadratic form in F,G,R and H. All 
scalar products can be expressed in the 5 independent variables s^, sj, 0„., d\ and (ft, such that 

]T | T | 2 = 2G 2 F | V us | 2 M^ 2 J 5 (^, s h 9„, d h <f>) . (5.14) 

spins 
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Carrying out the integrations over the remaining 4-3 — 5 = 7 variables in ( 5.13|) gives [|2 
dT 5 = G 2 F | V us | 2 N(s n ,si)J 5 (s n ,si,6 7T ,6i,(f))ds 7T dsid(cos8 7T )d(cos6i)d4> (5.15) 

where 



N{s w , 8l ) = (1 - Zl )a n X/(2 u 7i b M b K ) . 



(5.16) 



The form factors F, G, R and H are independent of and It is therefore possible to carry 
out two more integrations in (|5.15|) with the result 



dT 3 = G 2 F | V us | 2 N(s 7T ,s l )J 3 (s 7T ,s h 8 7T )ds 7T dsid(cos8 7T ) 
The explicit form of J5 is 



(5.17) 



■h 



(PL) 2 - (PN) 2 - s n s t + m 2 s n 



+ \G\ 2 (QL) 2 -(QN) 2 -Q 2 Sl + m 2 Q 2 



+ \R\ mj 
1 



+ 



M' K 



s t - m l 

(m 2 - Si) [Q 2 X 2 + s n (QL) 2 \ - < LNPQ > 2 

+ (F*G + FG*) [(PL) (QL) - (PN) (QN)} 
+ (F*R + FR*) m 2 [(PL) - (PN)} 

1 



K 



(F*H + FH*) (QN)(PL) 2 - (QL)(PL)(PN) - s nSl (QN) + mfs n (QL) 



+ (G*R + GR*) m] [(QL) - (QN)} 
1 



—y(G*H + GH*) (PL)(QL)(QN) - (PN)(QL) 2 + Sl (PN)Q 2 -mf(PL)Q' 
M K 

i 

M 



+ 



L K 



< LNPQ > -(F*G - FG*)M 2 K + (F*H - FH*)(PN) 



+ (G*H - GH*)(QN) + (R*H - RH*)m 2 



(5.18) 



For data analysis it is useful to represent this result in a still different form which displays 



the 9i and dependence more clearly |43| : 

J 5 = 2(1 — z{) 1% + h cos 29i + I 3 sin 2 81 ■ cos 20 + 7 4 sin 29i ■ cos + ^5 sin 81 • cos ( 
+ Iq cos 81 + ^7 sin 8[ ■ sin + J 8 sin 26*^ ■ sin + J 9 sin 2 ■ sin 20 . 

One obtains 

h = i{(l + ^)l^i| 2 + ^(3 + ^)(|F 2 | 2 +|F3| 2 )sin 2 ^ + 2^|F4| 2 } 



(5.19) 



-\(l-Zi){[F 1? -\ 



F 2 | 2 + |F 3 | 2 sin 2 ^ 
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h 
h 
h 
h 

h 



i(i-^){i^r-i^i 2 } 

z x ) Re(F*F 2 )sin# 



sin 2 6^ 



(I- 



{ Re(F*F 3 ) + z x Re(F 4 *F 2 )}sin^ 
{ Re(F 2 *F 3 ) sin 2 B v - z t Re(F*F 4 )} 
{ lm(F*F 2 )+ Zl Im(F 4 *F 3 )}sin^ 



■z t ) Im(F*F 3 )sin^ 
h = -^(l-^)M^3)sin 2 ft 



where 



F 3 
F 4 



7T ; 



X-F + a^PL) cos^-G 



(5.20) 



Mi 



■(PL)F -siR-a n X cosO^-G 



(5.21) 



The definition of Fi,...,F 4 in ( |5.21| ) corresponds to the combinations used by Pais and 
Treiman ||3| (the different sign in the terms ~ PL is due to our use of the metric diag(H — 

)). The form factors I\,...,Ig agree with the expressions given in ||43|| . We conclude 

that our convention for the relative phase in the definition of the form factors in Eq. ( |5.12|) 
agrees with the one used by Pais and Treiman. The comparison of ( |5.18 ) with j|4], table II] 
shows furthermore that it also agrees with this reference. 

The quantity J 3 can now easily be obtained from ( 5.19Q by integrating over <fr and 9i, 



J 3 = J dcj) d(cos9i)J5 = 8n(l — zi) 



h 



(5.22) 



5.4 Isospin decomposition 



The K14 decays ( |5.2j ) and ( p.3| ) involve the same form factors as displayed in Eq. ( |5.12| ) . We 



denote by A 00 and A _ the current matrix elements of the processes (|5.1|) - (|5.3| )- These 
are related by isospin symmetry 

A+- = ^ - A 00 . (5.23) 

This relation also holds for the individual form factors, which may be decomposed into a 
symmetric and an antisymmetric part under t <-> u (pi <-> p 2 ). Because of Bose symmetry 

5 We use the Condon-Shortley phase conventions. 
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and of the AI = | rule of the relevant weak currents, one has 

(F, G, R, H)oo = —(F + ,G~,R + ,H ) + _ 

(F,G,R,H) _ = V2(F-,G + ,R-,H + ) + _ (5.24) 

where 

F$_ = X - [F(s n , t, u) ± F(s n , u, £)] (5.25) 

and F(s n ,t,u) is defined in Eq. ( [5.121 ). 

The isospin relation for the decay rates is 

T(K + -> TT+TT-l+Vt) = —T(Kl -> 7T 7T ± / T Z/) + 2T{K + -> ttW+h) . (5.26) 

Isospin violating contributions affect the matrix elements and phase space, as a result of 
which this relation is modified. In order to illustrate the (substantial) effects from asymme- 
tries in phase space, we take constant form factors F, G and set R = 0, H — 0. Eq. ( [5.26 ) 
then reads (with physical masses for K + — > 7r + 7r~/ + z//, 7r°7r°/ + z/; and with M„.o = M n ± = 137 
MeV in K L -> tiV^z/) 

(16.0F 2 + 3.1G 2 )r = (20.1F 2 + 2.0G 2 )r 

T = ^,.10 2 sec _1 (5.27) 

in the electron mode and 

(1.79F 2 + 0.25G 2 )r = (2.64F 2 + 0.20G 2 )r (5.28) 

in the muon mode. 



5.5 Partial wave expansion 

The form factors may be written in a partial wave expansion in the variable 8 n . We consider 
a definite isospin Tin state. Suppressing isospin indices, one has [47], fi"S|] 

a w PL 



F 
G 
R 
H 



where 



z=o 

oo 



^/(cos^* 



X 



P l ( cos 9n)ri+ ~T7~" cos 6^(3 

oo 

£P/(cos^ 
z=o 



■ cos 9 W G 



(5.29) 



(5.30) 
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Table 5.1: Rates of K e ± decays. 





Jj events 






branching 


Particle Data 


DAFNE 


improve- 




ratio 


Group 


1 yr 


ment 


K + — ► 7T + 7r e + u e 


3.91 ■ 10~ 5 


3- 10 4 


3- 10 5 


10 


K + _^ 7T o 7r o e+z/e 


2.1 • 10" 5 


< 50 


2 • 10 5 


> 4 • 10 3 


i^i — > 7T°7i ± e T u 


6.2 ■ 10" 5 


16 


7- 10 4 


4- 10 3 



The partial wave amplitudes fi,gi,ri and hi depend on s n and s/. Their phase coincides 
with the phase shifts 5( in elastic scattering (angular momentum /, isospin J). More 
precisely, the quantities 

e- lS ^X 2l+l ; I = 0, 1, . . . ; X = f, g, r, h (5.31) 

are real in the physical region of decay. The form factors F\ and F4 therefore have a 
simple expansion 

F 1 = X^^iCcos^)/, 
1 

F A = -J^Piicose^iPLft + stn). (5.32) 
1 

On the other hand, the phase of the projected amplitudes 

F a = J Pz(cos^)F 2 rf(cos^) ; I = 0,1,2,... (5.33) 

is not given by S( , e.g., e~ l5l F2o is not real in the isospin one case. A similar remark applies 
to F 3 . 



5.6 Previous experiments 



We display in table [57L] the number of events collected so far. The data are obviously 



dominated by the work of Rosselet et al. 



which measures the n + n final state with 



good statistics. The authors parametrize the form factors as 

F = f s e< + / p e i,5 i 1 cos^ + D-wave 



G 
H 



ge 1 + D-wave 



he 



i6] 



D-wave 



(5.34) 
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with f s , f p , g and h assumed to be real [] . Furthermore, they put m e = 0, such that the 
form factors R and F 4 drop out in the decay distribution. Despite the good statistics, the 
experiment has not been able to separate out the full kinematic behaviour of the matrix 
elements. Therefore certain approximations/assumptions had to be made. For example, no 
dependence on S[ was seen within the limits of the data, so that the results were quoted 
assuming that such a dependence is absent. Similarly, f p was found to be compatible with 
zero, and hence put equal to zero when the final result for g was derived. A dependence on 
s n was seen, and found to be compatible with 

f s (q 2 ) = / s (0)[l + A /g 2 ] 
g(q 2 ) = g(0)[l + X 9 q 2 } 
% 2 ) = h(0)[l + X h q 2 ] 

q 2 = ( S7T -AM 2 )/AM 2 (5-35) 

with 

\f = \ g = \h = A. (5.36) 

These approximations to the form factors do not agree completely with what is found in 
the theoretical predictions. Dependence on si and non-zero values for higher partial waves 
all occur in the theoretical results. 

The experimental results for the threshold values and the slopes of the form factors are 

M 



f 3 (0) = 5.59 ±0.14 

0(0) = 4.77 ±0.27 

h(0) = -2.68 ±0.68 

A = 0.08 ±0.02. (5.37) 

We have used U | V us | = 0.22 in transcribing these results. (We note that from Eqs. ( |5.34j 
- 15.37| ) and f p = we obtain T Ke4 = (2.94 ± 0.16) • 10 3 sec -1 . This value must be compared 
with Tx eA = (3.26 ± 0.15) ■ 10 3 sec -1 obtained in the same experiment.) In addition to 
the threshold values ( [5.37D of the form factors, the phase shift difference 5 = 5$ — b\ was 
determined in five energy bins. The S-wave scattering length a[j was then extracted by 



using a model of Basdevant, Froggatt and Petersen []49 |. This model is based on solutions 



to Roy equations. The result for the scattering length is 

a° = 0.28 ±0.05. (5.38) 

A study by jp0"| , based on a more recent solution to Roy equations, gives 

a° = 0.26 ±0.05. (5.39) 
6 Note that, according to what is said in the previous subsection, the terms denoted by "D-wave" in Eq. 



(5.34) all contain (complex) contributions which are proportional to ^(cos^),/ > 0. 
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Turning now to the other channels, we consider the measured branching ratios 



BR(K + -> 7r°7T e + z/ e ) 



2 .o;°i 



• 10" 5 [25 events] [37] 



(5.40) 



and 



BR(K L -> n\ T e ± u) 



(6.2 ±2.0) • 10~ 5 

(5.8 ±0.2 ±0.4) ■ 10~ 5 



[16 events] [51] 
[780 ± 40 events] H53 



(5.41) 



The kinematic dependence of the form factors on the variables s w , si and 9 n has not yet been 
resolved experimentally in these decays. In order to proceed, we assume that the A 00 and 
v4 - form factors are independent of 9 n , e.g., F o = F 00 (s n ,t ± u) etc. As a result of this 
assumption, Goo, Hqq, Fo- and Rq- all vanish by Bose statistics. The contribution from i? o 
is completely negligible in the electron mode, and the contribution from the anomaly form 
factor to the decay ( |5.41|) is tiny. We neglect it altogether, as a result of which the above 
decays are fully determined by F o and Go-- We write 



oo 



and obtain for the rate 



2T 



K+- 



K L - 



>7r o 7r o e + I/e 



F (l±Ag 2 ), G - = G (l±Ag 2 ) 



|FoK,| 2 (2.01 ± 1.7A ± 0(A 2 )) • lO^ec" 1 
|G Ks| 2 (0.406 ± 0.47A ± 0(A 2 )) ■ lO^ec -1 



(5.42) 



(5.43) 
(5.44) 



where we have used physical phase space in flyp and Ml = M± = 137 MeV in §Aj ). 
This finally gives with A = 0.08 from ( 5.37|) 



l^ol 
\G Q \ 



5.5818;; [37] 

f 7.5 ±1.2 
I 7.3 ±0.3 





HI 



(5.45) 



which compares rather well with the isospin predictions (|5.24j) 



l^ol 
I Co | 



l/-(0)| 
y/2\g 



5.59 ±0.14 
= 6.75 ±0.38 



(5.46) 



5.7 Theory 



The theoretical predictions of K14 form factors have a long history which started in the sixties 
with the current algebra evaluation of F, G, R and H. For an early review of the subject 
and for references to work prior to CHPT we refer the reader to |26j (see also [27j). Here we 
concentrate on the evaluation of the form factors in the framework of CHPT |45|, [46|. We 
restrict our consideration to the isospin symmetry limit m u = ma, a = 0. 
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A) The one-loop result 



In Ref. [45, 46|, the form factors F, G and H have been evaluated in CHPT at order p 4 . 



The analytic expression for R has not yet been worked out to this accuracy |]5"3" |. This form 
factor only contributes significantly to decays. 

The chiral representation of the form factors at order E 2 was originally given by Weinberg 



a, 



F 
H 



G 



0. 



M. 



K 



V2F n 



3.74 



(5.47) 



We write the result for F at next-to-leading order in the form 

F(s n ,t,u) = ^f r {l + F + (s^,t, M ) + F-(^,t, 
F ± (s 7r ,t,u) = U^{s 7T ,t,u) + P^{s 7T ,t,u) + C^ 



u 



0(£ 4 )} 



(5.46 



and will use below an analogous expression for the form factor G. The superscript +(— ) 
denotes a term which is even (odd) under crossing t «-> u. The contributions Up(s n ,t,u) 
denote the unitarity corrections generated by the one-loop graphs which appear at order E A . 
They have the form 



Up(s n ,t,u) 



F- 2 [AqK) + a F (t) + a F (u)} 
F- 2 [b F (t) - b F (u)} 



(5.49) 



with 



A (s 7r ) 
a F (t) 



~Jkk( S tt) + „ 7r ^nr)( S 7r) 



4 



+ 



32 
1 

16 
1 



14MI + um 2 - m)r Kn (t) + ml + 2M 2 - 3t)r vK (t) 



{3M 2 K - 7M 2 + bt)K Kv (t) + {M 2 K - 5M 2 + 3t)K„ K (t) 



b F {t) 



9(L K7T (t) + L vK (t)) + (3M 2 K - 3M 2 - M)(M r Kw (t) + M r vK (t)) 



a F (t)--(M 2 K + M 2 -t)J r K7T (t). 



(5.50) 



The loop integrals ./^(s^), . . . which occur in these expressions are listed in the appendix 
|B[ The functions J P q and M P q depend on the scale \x at which the loops are renormalized. 
The scale drops out in the expression for the full amplitude (see below). 

The imaginary part of F 7 ^ 2 A (s 7r ) contains the / = 0, >S-wave tttt phase shift 



S (s n ) = (327rF 2 )-\2s n - M 2 )a n + 0(E 4 



(5.51) 
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as well as contributions from KK and r]f] intermediate states. The functions ap(t) and bp{t) 
are real in the physical region. 

The contribution P^ : (s 7r , t, u) is a polynomial in s n , t, u obtained from the tree graphs at 
order E A . We find 

1 9 

pp(s n ,t,u) = -^J2pfA s ^^^ u ) Lr i ( 5 - 52 ) 



F 2 

1 IT 1=1 



where 



32(^ - 2Ml) 
8(M| + s w - s,) 
2{M 2 K - %Ml + 5^ - si) 
32M^ 
AM 2 
2s t 

-2(t-u). (5.53) 

The remaining coefficients pf F are zero. The symbols L[ denote the renormalized coupling 
constants discussed in chapter |l[ 

Finally we come to the contributions Cp which contain logarithmic terms, independent 
of s n , t and u: 



Pl,F 
P2,F 

pIf 
pXf 
Pt,F 

P9,F 

PIf 



C 



F = (256ttX 2 



2 v2\-l 



hMi In ■ 



Ml 



Ml 



Ml 



2Ml ln^ -3ikL 2 in^ 



fi 2 



Cp = . 

The corresponding decomposition of the form factor G, 

G ± = U± + P± + C±, 

has the following explicit form: 

U£(s v ,t,u) = F- 2 [Ai(s^) + a G (t) + a G (u) 



U G (s^,t,u) 



K 2 Mt) - b G (u) 



(5.54) 
(5.55) 

(5.56) 



with 



a G (t) 



+ 



32 
1 

16 

3 



(2M| + 2M 2 n + 3t)J r Kn (t) - (2M| + 2M 2 - 3t) J^(t) 



-3M£ + 7M 2 - 5t)K Kv (t) + {-M 2 K + 5M 2 - 3t)i^(t) 



W*) + L vK (t) - (M 2 K - Ml + t)(M r K7T (t) + M r vK (t)) 
a G (t)-~(M 2 K + Ml-t)r K7r (t). 



(5.57) 
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The imaginary part of F n 2 Ai(s 7r ) contains the 1=1, P-wave phase shift 

5l(s n ) = (dGrrF^s^J 2 + 0(P 4 ) (5.58) 

as well as contributions from KK intermediate states. The functions cig, be are real in the 
physical region. 



The polynomials 



^4Ei(^M)i: (5.59) 

7T i=l 



arc 



p+G = ~2(M 2 K + s n - Sl ) 

Ptc = 4M* 

Ptc = 2s i 

P%g = S(t-u) 

P3,G = \P2,G- ( 5 - 6 °) 

The remaining pf G vanish. The logarithms contained in Cq are 

C% = -C$. (5.61) 
The form factor H starts only at 0(E 4 ). The prediction is 

# = -^t = -2.66 (5.62) 

in excellent agreement with the experimental value. 

The results for F and G must satisfy two nontrivial constraints: i) Unitarity requires 
that F and G contain, in the physical region 4M 2 < s n < (Mk — rrii) 2 , imaginary parts 
governed by S- and P-wave nir scattering [these imaginary parts are contained in the func- 
tions A (s 7r ), A 1 (s 7r )]. ii) The scale dependence of the low-energy constants L\ must be 
compensated by the scale dependence of Up,G an d CV,g f° r au values of s^, t, u, M 2 , M\. 
[Since we work at order E A , the meson masses appearing in the above expressions satisfy the 
Gell-Mann-Okubo mass formula.! We have checked that these constraints are satisfied. 



B) Comparison with experiment 



One striking feature of the chiral prediction for the form factors is that the only important 
dependence on the low-energy constants is through L\, L 2 and L 3 . We proceed by fixing 
L4, L5 and Lg at the values found in other processes (see chapter J], in particular table 1). 

A rather extensive analysis of the chiral prediction and the data on K + — > 7r + 7r~e + i/ and 
elastic tttt scattering has been given in Refs. |45|, JRJ. We refer the reader to these articles 
for details. Here we mention the following points. 
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Table 5.2: Predictions of chiral symmetry following from the fit to the K e ± data [H alone 
(column 3) and the combined determination from txtx [55 and K ei data [Q (last column). 
The first column gives the prediction of the leading order term in the low-energy expansion 
of the txtx amplitude. 





leading 










order 


experiment 


K ei alone 


K ei ± TXTX 






0.08 ±0.02 


0.06 ±0.02 


0.06 ±0.02 




0.16 


0.26 ±0.05 


0.20 


0.20 


h° 


0.18 


0.25 ±0.03 


0.26 


0.26 


a 2 Q 


-0.045 


-0.028 ±0.012 


-0.040 


-0.041 


bl 


-0.089 


-0.082 ±0.008 


-0.069 


-0.070 


a\ 


0.030 


0.038 ±0.002 


0.037 


0.036 


b\ 






0.045 


0.043 


al 




(17 ±3) • 10" 4 


21 • 10" 4 


20 • 10" 4 


a 2 2 




(1.3 ±3) • 10" 4 


3.5 • 10~ 4 


3.5 • 10~ 4 



1. Fixing Li, L2 and L3 from / s (0), g(0) and the slope A/ of the form factor f s gives 

L[(M p ) = [0.5 ±0.3] • 10~ 3 
L r 2 (M p ) = [1.6 ±0.3] • 10~ 3 

L 3 = [-3.2 ±1.1] -10~ 3 . (5.63) 



The error bar corresponds |46|] to an increase of x 2 by one. It does not include the error 
due to unknown higher order corrections in the chiral expansion of the form factor. 
Having determined Li,L 2 and L 3 , one may then work out the form factors from the 
representation ( |5.48| - |5T6T|) . The result is shown in Fig. |5.2| , where we plot for the 
electron mode the quantity 



1 



dsi 



d(cos9 n )F(s n , si, cos9 n ) 



2>> 1/2 



(5.64) 



and similarly for g{s. K ). The lowest order results Eq. ( |5.47p (labelled "tree") plus the 
experimental central values and the central values corresponding to Eq. ( |5.63|) are 
displayed. Note that the slope of the g form factor has not been included in the fit 
and is thus a prediction. It matches very well with the experimental data. 



2. The decay K 



7r + 7r e + u e allows one to test the large Nc prediction 
(L r 2 -2Ll)/L 3 = (large N c ). 



(5.65) 



From the values in Eq. ( |5.63| ), we see that a small non-zero result for this combination 
is preferred, but that it is consistent with zero within the errors. The fit was also done 
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Figure 5.2: The form factors / s (s,r) and g(s 7r ) (Eq. |5.64|) according to the chiral representa- 
tion (electron mode). The dotted lines show the lowest order result ( |5.47|) , and the dashed 
lines correspond to L ly L 2 and L 3 from ( |5.63|) . The experimental result ( |5.37|) is displayed 
by a solid line. 



46 1 using the variables 



X\ — L 2 — 2_Z^ — L 3 
X2 = L 2 

X 3 = (L r 2 -2Ll)/L 3 (5.66) 



with the result 



X x = (3.8 ±0.9) • 1(T 3 
X 2 = (1.6 ±0.3) • 1(T 3 

X 3 = -0.19+°if . (5.67) 

(Xi and X 3 are scale independent, X 2 is evaluated at the rho mass.) The result is that 
the large Nc prediction works remarkably well. 

3. Having determined the low-energy constants, one is in a position to study the predic- 
tions. The coefficients Lx,L 2 and L 3 also govern elastic 7T7r-scattering, and the real 
test of the theory is that these coefficients are simultaneously compatible with the 
elastic tctt amplitude. The most straightforward way to check this is to predict the nir 
threshold parameters. The chiral predictions were worked out in Ref. |5(| |57f . If we 
use the determination ( |5.63|) , we obtain the prediction in table |5]^, third column. (For 



l 3 , £4 which occur in af , b\ we have used the central value l 3 = 2.9, Z 4 = 4.3 from Ref. 
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6|). The predictions are within l| standard deviations of the data in all cases. Note, 



in particular, the nice agreement for the I = 0, 2 D-wave scattering lengths a 2 , 



Furthermore, it is comforting to see that the SU(2) xSU(2) prediction []5B|, pTF 



0.20 ±0.01 



(5.66 



survives the K e4 test unharmed. 



4. It is of interest to provide the best determination of the low-energy constants by in- 
cluding the maximum amount of data. This includes the K e 4 form factors / s (0),^(0) 
and Xf, as well as the direct determination of 5® — S\ in K ei decay. We take the other 
information as the tttt threshold parameters well as the universal curve 



a l) 



2a° - 5a 2 - 0.96(a£ 
-0.7(a°-0.3) 2 
0.69 ±0.04 . 



0.3) 



(5.69) 



The results of the fit are shown in the last column of table f>.2[ The corresponding 
values for L\, L 2 and L 3 are 



L\{M P ) 
L r 2 (M p ) 



(0.7 ±0.5) • 10~ 3 
(1.2 ±0.4) • 10" 3 
(-3.6 ±1.3) • 10~ 3 



(5.70) 



The error includes the theoretical error bar, see Ref. [pq| . [The one-loop representation 
(|5.48 - |5.62 ) of the form factors F, G and H, evaluated at the central values ( |5.70|) , gives 
r^ e4 = 2.5 • 10 3 sec -1 . This is somewhat lower than the experimental 0] width 



3.26 



13 

10 3 



sec 1 and the value Y 



2.94 • 10 3 sec 1 which follows from 



( |5.37|) . The reason for considering nevertheless ( |5.70| ) as the present best estimate for 
Li,L 2 and L 3 is discussed at some length in |f4B 
tree level result F 



since the chiral corrections to the 



G = 3.74 are large, one should not expect that the one- loop 
corrections already do the complete job - rather, higher order terms have any right to 
also contribute accordingly. The above result for Li,L 2 and L3 includes an estimate 
for these additional terms. (The experimental width for K ei is within the uncertainties 
ALi, AL 2 and AL 3 quoted.)] 



5.8 Improvements at DAFNE 

The chiral analysis of decays has been used so far for three purposes: 



1. The K e4 data from Ref. [[Pj make predictions for the slope of the G form factor and 
for the 7T7T scattering lengths. These are given in table |5.2| . 
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2. The same K e4: data allow one to test the large Nc prediction, see Eqs. ( |5.65| - |5.67| ). 



3. The full set of K e4 and irir scattering data allows the best determination of the coeffi- 
cients Li,L 2 and L 3 in the chiral Lagrangian, see Q5.7C ). 



In the next generation of A^ 4 decay experiments, there is the opportunity to improve the 
phenomenology of K i4: (see table |5.1| ): 



1. The present experimental uncertainty on G is still too large to provide a precise value 
for the large Nq parameter (L 2 — 2L r l ) / L3. (K° — > 7r°7r _ e + z/ e decays are mainly sensitive 
to G\_ which in turn can be used to pin down L3. K + — > 7r°7r°e + z/ e is mainly sensitive 
to F+_ which contains L±,L 2 and L3.) 

2. The observation of all K14 reactions with high statistics could provide a cleaner sepa- 
ration of the various isospin amplitudes. 

3. A very useful innovation would be to analyze the experimental data directly using the 
framework of chiral perturbation theory Rather than making assumptions about the 
absence of P- waves, D-waves etc., one could parametrize the data using the full chiral 
perturbation formulas, and directly decide the quality of the fit and the favoured values 
of the low-energy constants. 

4. Finally, we come to a most important point. As we mentioned already, K + — > 
7r + 7r~e + i/ e has been used |50| to determine the isoscalar S'-wave scattering length with 



the result a{] = 0.26 ± 0.05. This value must be compared with the SU(2) xSU(2) 
prediction [|56], [57]] a[j = 0.20 ± 0.01. Low-energy tttt scattering is one of the few places 
where chiral symmetry allows one to make a precise prediction within the framework 
of QCD. In their article, Rosselet et al. comment about the discrepancy between 
a[] = 0.26 ± 0.05 and the leading order result ]59| a[j = 0.16 in the following manner: 



"... it appears that this prediction can be revised without any fundamental change in 
current algebra or in the partial conservation of axial- vector current |K], |6H|." Today, 
we know that this is not the case: It would be a major difficulty for QCD, should the 
central value a° = 0.26 be confirmed with a substantially smaller error. 

decays are - at present |62|] - the only available source of clean information on tttt 



S'-wave scattering near threshold. We therefore feel that it would be rather appropriate 
to clarify this issue. 
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A Notation 



The notation for phase space is the one without the factors of 2ir. For the decay rate of a 
particle with four momentum p into n particles with momenta p\, . . . ,p n this is 

n n d 3 pi 

dups(p;pi, ■ ■ -,Pn) = 8 4 {p - ^Pi) n iri ■ ( AA ) 

We use a covariant normalization of one-particle states, 

<p'\p>= (27i) 3 2p°5 3 (p' -p) , (A.2) 

together with the spinor normalization 

u(p,r)u(p, s) = 2m5 rs . (A. 3) 

The kinematical function X(x,y, z) is defined as 

X(x, y, z) = x 2 + y 2 + z 2 — 2{xy + yz + zx) . (A. 4) 

We take the standard model in the current x current form, i.e., we neglect the momentum 
dependence of the VT-propagator. The currents used in the text are : 



1 



Q = diag(2/3, -1/3, -1/3) . (A.5) 

The numerical values used in the programs are the physical masses for the particles as given 
by the Particle Data Group QX[ . In addition we have used the values for the decay constants 
derived from the most recent measured charged pion and kaon semileptonic decay rates Jl], [18| 

= 93.2 MeV 

F K = 113.6 MeV. (A.6) 

We do not need values for the quark masses. For the processes considered in this report we 
can always use the lowest order relations to rewrite them in terms of the pseudoscalar meson 
masses (see chapter [I]). For the KM matrix element |V^ S | we used the central value, 0.220, 
of Ref. jl]. The numerical values for the L[(M p ) are those given in chapter [I]. 

The number of events quoted for DAFNE are based on a luminosity of 5 • 10 32 cm _2 s _1 , 
which is equivalent [Q to an annual rate of 9 ■ 10 9 (1.1 • 10 9 ) tagged K ± {K£) (1 year = 
10 7 s assumed). 

Whenever we quote a branching ratio for a semileptonic K° decay, it stands for the 
branching ratio of the corresponding Kl decay, e.g., 

BR(K° -> ir~l + v) = BR(K L -> tt ± / t z/) . (A.7) 

We use the Condon-Shortley phase conventions throughout. 
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B Loop integrals 



In this appendix we define the functions appearing in the loop integrals used in the text. 
First we define the functions needed for loops with two propagators, mainly in the form 
given in Ref. [35|. We consider a loop with two masses, M and m. All needed functions can 
be given in terms of the subtracted scalar integral J(t) = J(t) — J(0), 

f d 4 p 1 
J ^ = V (2ir) d ((p + kf - M 2 )(p 2 - m 2 ) ^'^ 

with t = k 2 . The functions used in the text are then : 

, 1 f 1 , , M 2 — txil -x)-Ax 



16tt 2 7 M 2 — Ax 

1 f A m 2 E, m 2 , (t + vT) 2 - A 2 

' 2 H log — - — — log — log ■ 



32vr 2 [ t M 2 A M 2 t to (t - v^)2 - A 2 j ' 
J r (t) = J (t) -2k, 

M r (t) = —{t-2X}J(t)+=^J(t) + —^-- 
w 12* 1 1 W 3t 2 w 288tt 2 6 

1 f M 2 m 2 , m 2 ] 
-96rt{ E + 2 — l0g Mj| ■ 

A 2 - 

L(t) = —J(t), 
K{t) = pit), 

H(t) = 2 -^ 2 t + jr 2 [tM r (t)-L(t)], 

A = M 2 - m 2 , 
E = M 2 + m 2 , 

A = A(t,M 2 ,m 2 ) = (t + A) 2 -4*M 2 . (B.2) 

In the text these are used with subscripts, 

Jy(i) = J(t) with M = Mi,m = Mj (B.3) 

and similarly for the other symbols. The subtraction point dependent part is contained in 
the constant k 

1 M 2 log(^)-m 2 log(^) 



k = — v M 7 VM 7 fB 4) 

32vr 2 M 2 -m 2 ' 1 J 

where /x is the subtraction scale. 

In addition, in section |] these functions and symbols appear in a summation over loops 
/ with 

J T (t) = J(f) with M = M I ,m = m I ; 

£/ = Mf + m] (B.5) 
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and again similarly for the others. There the combination B 2 appears as well : 

B 2 (t,M 2 ,m 2 ) = B 2 (t,m 2 ,M 2 ) (B.6) 

1 t) X(t,M 2 ,m 2 )J(t) | tS-8M 2 m 2 lQg M 2 



288tt 2 v ' 12t 384tt 2 A ° m 2 

The last formula to be defined is the three propagator loop integral function C(ti,t2, M 2 , m 2 ) 
where one of the three external momenta has zero mass and two of the propagators have the 
same mass M. Here t± = (qi + q 2 ) 2 , t 2 = q\ and q\ = 0. 



C(t u t 2 ,M 2 ,m 2 ) = -if 



d 4 p 1 



(27r) d (p 2 - M 2 )((p + qi ) 2 - M 2 )((p + qi + q 2 ) 2 - m 2 ) 

16^1 dX Jo dy M 2 - y(A + tT) + xyfa - t 2 ) + y 2 t x 
Li 2 ( — ttt I + Li 2 



(4 7 r) 2 (t 1 -t 2 ) [ '\y + (t 2 )J A \y-(t 2 ) 

-Li 2 ( — ]— ) - Li 2 ' ' 



V±(t) = ^{t + Aiy^M 2 ,™ 2 )} (B.7) 



where Li 2 is the dilogarithm 



Li 2 (x) = - C — log(l - xy) . (B.8) 
Jo y 
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C Decomposition of the hadronic tensors 1^ 

Here we consider the tensors 

Jf" = / dxe^+^v < | TV» m (x)I»_ l5 (y) \ K + (p) > , I = V,A (C.l) 



and detail its connection with the matrix element (1.2). 

The general decomposition of A^, V^ v in terms of Lorentz invariant amplitudes reads 
[g |] for q V 

+ Ai(qWg^ u - W^q u ) + A 2 {q 2 g'" u - q^q u ) 
\2F K {F${q 2 )-l 



and 



(Ml - H/ 2 )g 
1 



+ A 3 } (qWq^ - q 2 W^)W u (C.2) 



^=V^ = iV^q^p (C.3) 

where the form factors Ai(q 2 ,W 2 ) and Vi(q 2 ,W 2 ) are analytic functions of q 2 and IV 2 . 
Fy(q 2 ) denotes the electromagnetic form factor of the kaon (Fy(0) = 1). A^ v satisfies the 
Ward identity 

q^ = -V2F K p». (C.4) 

In the process ( [Lip the photon is real. As a consequence of this, only the two form 
factors A^O^ 2 ) and Vi{Q,W 2 ) contribute. We set 

A(W 2 ) = Ax(0,W 2 ) 

V(W 2 ) = Vx(0,W 2 ) (C.5) 



and obtain for the matrix element ( |1.2| ) 

T = -iG F /V2eV u *el {V2F K l? - (V" - A^)lJ] , (C.6) 



9^=0 



with 



K = F + m l u(p v )(l + l5 ) ^l±^— v( Pl ) . (C.7) 

m i - ipi + Q) 



Grouping terms into an IB and a SD piece gives ( |0| , |L~B| ) . As a consequence of ( |C4| ), T is 
invariant under the gauge transformation — > e M + g M . 

The amplitudes t4 1; v4 2 and Vi are related to the corresponding quantities Fa, R and Fy 
used by the PDG [|TJ by 

- y/2M K (A h A 2 , K) = (F A , R, F v ). (C.8) 
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The last term in ( |C.2| ) is omitted in PJ. It contributes to processes with a virtual photon, 

Finally, the relation to the notation used in |3| is 

2(A±V) 2 = (a k ±v k ) 2 [2] 

V2(A,V) = (Fa, Fy) § . (C.9) 
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D Formulas for the traces in terms of x, y and z for the 
decays K+ -> l + ul' + l'- 

This is the FORTRAN program used to evaluate the differential decay rate in terms of 
the kinematic variables used in the text. The formfactors are Al,A2,A4 and VI with their 
complex conjugates A1C, A2C, A4C and VIC, all made dimensionless by multiplying with 
the relevant power of Mk- T is the quantity j— ^2 spins T M *T M 1. The matrix element squared, 
before the integration over the lepton pair kinematic variables, is available on request from 
the authors. 

C ALL QUANTITIES ARE IN UNITS OF MK TO THE RELEVANT POWER 
W2 = 1.0 - X + Z 
P$W = 1.0 - X/2.0 
P$PN = P$W - Y/2.0 
PL$PN= (W2 - RD/2.0 
PL$W = RL + PL$PN 
PN$W = W2 - PL$W 
Q$W = X/2.0 - Z 
Q$PL = Y/2.0 - PL$W 
Q$PN = P$PN - PN$W 
DEN0M1 = 1.0/(2.0*Q$PL + Z) 
DEN0M2 = 1.0/ (X - Z) 
All = REAL ( Al + A1C) 
A22 = REAL ( A2 + A2C) 
A44 = REAL ( A4 + A4C) 
VV = REAL ( VI + VIC) 
A12 = REAL (A1*A2C + A10A2) 
A14 = REAL (A1*A4C + A10A4) 
A1V = REAL (A1*V1C + AlOVl) 
A24 = REAL (A2*A4C + A20A4) 
A2V = REAL (A2*V1C + A20V1) 
A4V = REAL (A4*V1C + A40V1) 
A1A1C = REAL (A1*A1C) 
A2A2C = REAL (A2*A2C) 
A4A4C = REAL (A4*A4C) 
V1V1C = REAL (V1*V1C) 
T = 0.0 

T = T + A1A1C * ( 16*Q$PL*Q$PN*W2 - 16*Q$PL*Q$W*PN$W - 16*Q$PN* 
+ Q$W*PL$W - 8*PL$PN*Z*W2 ) 

T = T + A2A2C * ( - 16*Q$PL*Q$PN*Z - 8*PL$PN*Z**2 ) 

T = T + A4A4C * ( 8*Q$W**2*PL$PN*Z*W2 " 16*Q$W**2*PL$W*PN$W*Z - 
+ 8*PL$PN*Z**2*W2**2 + 16*PL$W*PN$W*Z**2*W2 ) 

T = T + V1V1C * ( - 8*P$PN*Q$PL*X + 8*P$PN*Y*Z + 16*Q$PL*Q$PN 
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+ - 4*Q$PN*X*Y ) 

T = T + FK**2*DEN0M1**2 * ( " 32*Q$PL**2*PL$PN*Z** (-1) *RL - 32* 
+ Q$PL*Q$PN*RL - 32*Q$PL*PL$PN*RL + 32*Q$PN*RL**2 + 8*PL$PN*Z* 
+ RL + 32*PL$PN*RL**2 ) 

T = T + FK**2*DEN0M1*DEN0M2 * ( 32*P$PN*Q$PL*RL + 32*Q$PL*PL$PN*X 
+ *Z**(-1)*RL - 16*Q$PN*Y*RL - 32*PL$PN*Y*RL ) 

T = T + FK**2*DEN0M2**2 * ( - 8*PL$PN*X**2*Z** (-1) *RL + 32*PL$PN 
+ *RL ) 

T = T + FK*DEN0M1*A11 * ( 16*Q$PN*Q$W*RL - 16*Q$PN*PL$W*RL + 16* 

+ Q$W*PL$PN*RL + 8*PN$W*Z*RL ) 
T = T + FK*DEN0M1*A22 * ( - 16*Q$PL*Q$PN*RL + 24*Q$PN*Z*RL + 16* 

+ PL$PN*Z*RL ) 
T = T + FK*DEN0M1*A44 * ( 16*Q$PL*Q$W*PN$W*RL - 8*Q$PN*Z*W2*RL + 

+ 8*Q$W*PN$W*Z*RL - 16*PL$W*PN$W*Z*RL ) 
T = T + FK*DEN0M1*VV * ( 16*P$PN*Z*RL - 8*Q$PN*X*RL ) 
T = T + FK*DEN0M2*A11 * ( - 16*P$PN*Q$W*RL + 16*P$W*Q$PN*RL ) 
T = T + FK*DEN0M2*A22 * ( - 16*P$PN*Z*RL + 8*Q$PN*X*RL ) 
T = T + FK*DEN0M2*A44 * ( 16*P$W*PN$W*Z*RL - 8*Q$W*PN$W*X*RL ) 
T = T + A12 * ( - 8*Q$PL*PN$W*Z - 8*Q$PN*PL$W*Z - 8*Q$W*PL$PN*Z 

+ ) 

T = T + A14 * ( 8*Q$PL*PN$W*Z*W2 + 8*Q$PN*PL$W*Z*W2 - 16*Q$W*PL$W 
+ *PN$W*Z ) 

T = T + A1V * ( - 8*P$PN*Q$PL*Q$W - 8*P$PN*PL$W*Z - 4*Q$PL*PN$W* 
+ X + 4*Q$PN*Q$W*Y + 4*Q$PN*PL$W*X + 4*PN$W*Y*Z ) 

T = T + A24 * ( 8*Q$PL*Q$W*PN$W*Z + 8*Q$PN*Q$W*PL$W*Z - 8*Q$W**2* 
+ PL$PN*Z + 8*PL$PN*Z**2*W2 - 16*PL$W*PN$W*Z**2 ) 

T = T + A2V * ( - 16*P$PN*Q$PL*Z + 8*Q$PN*Y*Z ) 

T = T + A4V * ( 8*P$PN*Q$PL*Z*W2 - 8*P$PN*Q$W*PL$W*Z - 8*P$W*Q$PL 
+ *PN$W*Z + 8*P$W*Q$PN*PL$W*Z - 4*Q$PN*Y*Z*W2 + 4*Q$W*PN$W*Y*Z 
+ ) 

T = -T 
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E FORTRAN routine for the calculation of the re- 
duced square of the i^/3 7 matrix element SM 



FUNCTION SM(EG,EP,W2,EL,X) 
C THE FUNCTION SM CALCULATES THE REDUCED SQUARE OF THE MATRIX 
C ELEMENT OF SECT. 4 IN TERMS OF THE SCALAR VARIABLES EG (PHOTON 
C ENERGY), EP(PI0N ENERGY), W2 (INVARIANT MASS SQUARED OF LEPTON 
C PAIR), EL (ENERGY OF CHARGED LEPTON) AND X=PL.Q/MK~2 AND IN 
C TERMS OF THE VECTOR AMPLITUDES B1,..,B7, AXIAL AMPLITUDES 
C Al, A2, A3 (A4=0 TO 0(P~4)) AND CI, C2 DEFINED IN SECT. 4. 
C ALL DIMENSIONFUL QUANTITIES ARE NORMALIZED TO THE KAON MASS: 
C MP=M(PI0N)/M(KA0N), ML=M (LEPTON) /M (KAON) , ETC. 

REAL ML,MP,ML2,MP2,ML4 

COMMON/MASSES/ML, MP 

ML2=ML**2 

MP2=MP**2 

ML4=ML**4 

C SCALAR PRODUCTS: QPP=Q.P', PLPP=P ( LEPTON ) .P' , WPP=W.P', 
C QW=Q.W, ALL SCALED TO M(KA0N)=1; W=P (LEPTON) +P (NEUTRINO) 

QPP=EP+EG+ (W2-MP2-1 . ) /2 . 

PLPP=EL-X- (W2+ML2) /2 . 

WPP=-EG+ ( 1 . -MP2-W2) /2 . 

QW=-EP+ ( 1 . +MP2-W2) /2 . 
C FOR ILLUSTRATION, THE TREE LEVEL AMPLITUDES B1,...,B7, 
C A1,A2,A3,C1,C2 FOR K0(L3GAMMA) ARE LISTED BELOW 

Bl=-1. 

B2=0. 

B3=2./QPP 

B4=0. 

B5=0. 

B6=1./QPP 

B7=B3 

A1=0. 

A2=0. 

A3=0. 

Cl=2. 

C2=l. 

C IN THE FOLLOWING, SM IS CALCULATED IN TERMS OF SCALAR 
C PRODUCTS, MASSES AND INVARIANT AMPLITUDES. 
C THIS PART IS INDEPENDENT OF THE CHOICE OF SCALAR VARIABLES 
C TO SPECIFY THE KINEMATICS (P (LEPTON). Q IS DENOTED X) 
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R1=B1*(B1-B5*ML2)*(-ML2 + W2) 
R1=R1 + B1*B2*2*ML2*(X-QW) 

R11=B1*B3+B2*B4*W2+(B2*B7+B3*B4)*WPP+B3*B7*MP2 

R1=R1+R11*(4*X*PLPP-2*X*WPP-2*PLPP*QW- 
+ QPP*ML2 + QPP*W2) 

R1=R1 + B1*(B4+B6)*2*ML2*(PLPP-WPP) 

R1=R1 + B1*A1*(-4*X*WPP + 4*PLPP*QW) 

R1=R1 + B1*A2*(4*X*W2-2*QW*ML2-2*QW*W2) 

R1=R1 + B2*(B2*W2+2*B3*WPP)*2*X*(X-QW) 

R1=R1 + B2*(B5*W2+B6*WPP)*2*ML2*(X-QW) 

R12=B2*A1+B3*A2+A3*ML2* (A2-C2/2 . /X) 

R1=R1+R12*(-2*X*QW*WPP+2*X*QPP*W2+2*PLPP 
+ *QW**2-QW*qPP*ML2-QW*QPP*W2) 

R1=R1 + B3**2*2*MP2*X*(X-QW) 

R1=R1 + B3*(B5*WPP+B6*MP2)*2*ML2*(X-QW) 

R1=R1+(B4**2*W2+B7**2*MP2+2*B1*B7) * (-1 . /2 . *ML2*MP2+1 . /2 . *MP2*W2 
+ + 2*PLPP**2-2*PLPP*WPP) 

R1=R1+ (B4* (B5*W2+B6*WPP) +B7* (B5*WPP+B6*MP2) ) *2*ML2* (PLPP-WPP) 

R1=R1+B4*B7*WPP*(4*PLPP**2-4*PLPP*WPP- 
+ ML2*MP2+MP2*W2) 

R1=R1+ (B4*A1+B7*A2) * (-2*X*WPP**2+2*X*MP2*W2+2*PLPP 
+ *QW*WPP-2*PLPP*QPP*W2-QW*ML2*MP2-QW*MP2*W2 
+ + QPP*WPP*ML2 + QPP*WPP*W2) 

R1=R1+ (B5* (B5*W2/2+B6*WPP) +B6**2*MP2/2) *ML2* (ML2-W2) 

R1=R1 + A1**2*(2*X**2*MP2-4*X*PLPP*QPP-2*X* 
+ QW*MP2 + 2*X*QPP*WPP + 2*PLPP*QW*QPP) 

R1=R1 + A1*A2*(-4*X**2*WPP + 4*X*PLPP*QW + 2*X 
+ *QW*WPP + 2*X*QPP*ML2-2*PLPP*QW**2-QW*QPP* 
+ ML2-QW*QPP*W2) 

R1=R1 + (Al+Cl/2./X)*A3*ML2*(2*X*QW*MP2-2*X*QPP*WPP 
+ -2*PLPP*QW*QPP-2*QW**2*MP2 + 4*QW 
+ *QPP*WPP + QPP**2*ML2-QPP**2*W2) 

R1=R1 + A2**2*(2*X**2*W2-2*X*QW*ML2-2*X*QW* 
+ W2 + QW**2*ML2 + QW**2*W2) 

R1=R1+A3**2*ML2* (1 . /2 . *QW**2*ML2*MP2-1 . /2 . *QW**2 
+ *MP2*W2-QW*QPP*WPP*ML2+QW*QPP*WPP*W2 +1./ 
+ 2 . *QPP**2*ML2*W2-1 . /2 . *QPP**2*W2**2) 

R2=C1**2*(-1./2.*ML4*MP2 + 1 . /2 . *ML2*MP2*W2 
+ + X**2*MP2-2*X*PLPP*QPP-X*QW*MP2 + 2* 
+ X*QPP*WPP-X*ML2*MP2 + 2*PLPP**2*ML2 + 2* 
+ PLPP*QPP*ML2-2*PLPP*WPP*ML2 + QW*ML2*MP2-2* 
+ QPP*WPP*ML2) 

R2=R2 + C1*C2*(2*X**2*WPP-2*X*PLPP*QW + 4*X* 



74 



+ PLPP*ML2-X*QPP*ML2 + X*QPP*W2-4*X*WPP* 
+ ML2 + 2*PLPP*ML4-2*WPP*ML4) 

R2=R2 + C2**2*(1./2.*ML**6-1./2.*ML4*W2 + 2*X**2* 
+ ML2 + X**2*W2-3*X*QW*ML2 + 2*X*ML4-X* 
+ ML2*W2-QW*ML4) 

RI=B1*C1*(-2*X*WPP + 2*PLPP*QW-2*PLPP*ML2 
+ + QPP*ML2-QPP*W2 + 2*WPP*ML2) 

RI=RI + B1*C2*(-ML4 + ML2*W2-2*X*ML2-2*X* 
+ W2 + 4*QW*ML2) 

RI=RI + B2*C1*(2*X**2*WPP-2*X*PLPP*QW-2*X* 

+ PLPP*ML2-2*X*PLPP*W2-2*X*QW*WPP-X*QPP* 

+ ML2 + X*QPP*W2 + X*WPP*ML2 + X*WPP*W2 + 2* 

+ plpp*qw**2 + plpp*qw*ml2 + PLPP*qw*w2 + qw*qpp* 

+ ML2-qW*qPP*W2 + l./2.*qPP*ML4-l./2.*qpp*W2**2) 
RI=RI+B2*C2*ML2*(-2*X**2-X*ML2-X*W2 

+ +2*qw**2+qw*ML2+qw*w2) 

RI=RI + B3*C1*(2*X**2*MP2-4*X*PLPP**2-4*X* 
+ PLPP*qpp + 2*X*PLPP*WPP-2*X*qW*MP2 + 2*X* 

+ qpp*wpp + 2*PLPP**2*qw + 2*PLPP*qw*qpp + PLPP*qpp* 

+ ML2-PLPP*qpp*W2 + qPP**2*ML2-qPP**2*W2) 

RI=RI + B3*C2*(-4*X**2*PLPP + 2*X**2*WPP + 2*X* 
+ PLPP*qW-2*X*PLPP*ML2-X*qpp*ML2-X*qPP* 
+ W2 + 2*PLPP*qW*ML2 + 2*qW*qpp*ML2) 

RI=RI + B4*C1*(1./2.*ML4*MP2-1./2.*MP2*W2**2 + 2* 
+ X*PLPP*WPP-2*X*WPP**2 + X*MP2*W2-2*PLPP**2 
+ *qw-2*PLPP**2*ML2-2*PLPP**2*W2 + 2*PLPP*qw*WPP- 
+ PLPP*qpp*ML2-PLPP*qPP*W2 + 2*PLPP*WPP*ML2 + 2* 
+ PLPP*WPP*W2-qW*MP2*W2 + qPP*WPP*ML2 + qPP*WPP*W2) 

RI=RI + B4*C2*ML2*(-2*X*PLPP+X*WPP- 
+ PLPP*qW-PLPP*ML2-PLPP*W2 + 2*qW*WPP 
+ +l./2.*qPP*ML2-l./2.*qPP*W2 + WPP*ML2 
+ + WPP*W2) 

RI=RI + B5*Cl*ML2*(X*WPP-PLPP*qW-PLPP* 
+ ML2-PLPP*W2-l./2.*qPP*ML2 + l./2.*qPP*W2 
+ +WPP*ML2 + WPP*W2) 

RI=RI + B5*C2*ML2*(-1./2.*ML4+1./2.*W2**2-X* 
+ ML2+qW*W2) 

RI=RI + B6*C1*ML2*(X*MP2-2*PLPP**2-2*PLPP 

+ *qpp+2*PLPP*wpp-qw*MP2 + 2*qpp*wpp) 

RI=RI + B6*C2*ML2*(-2*X*PLPP+X*WPP+ 
+ PLPP*qW-PLPP*ML2 + PLPP*W2-l./2.*qpp* 
+ ML2 + l./2.*qPP*W2) 

RI=RI + B7*Cl*(2*X*PLPP*MP2-X*WPP*MP2-4* 
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+ PLPP**3-4*PLPP**2*QPP + 4*PLPP**2*WPP-PLPP*QW* 
+ MP 2 + 4*PLPP*QPP*WPP + PLPP*ML2*MP2-PLPP*MP2* 
+ W2 + l./2.*QPP*ML2*MP2-l./2.*QPP*MP2*W2) 
RI=RI + B7*C2*(-4*X*PLPP**2 + 4*X*PLPP*WPP-X 
+ *MP2*W2-2*PLPP**2*ML2 + 2*PLPP*WPP*ML2 + QW*ML2 
+ *MP2) 

RI=RI + A1*C1*(2*X**2*MP2-4*X*PLPP*QPP + 2*X* 
+ PLPP*WPP-2*X*QW*MP2 + 4*X*QPP*WPP-X*ML2 

+ *MP2-X*MP2*W2-2*PLPP**2*QW + PLPP*QPP*ML2 + 
+ PLPP*QPP*W2 + 2*QW*ML2*MP2-2*QPP*WPP*ML2) 

RI=RI + A1*C2*(4*X**2*WPP-4*X*PLPP*QW-2*X* 
+ QPP*ML2 + 2*QW*QPP*ML2) 

RI=RI + A2*C1*(-2*X**2*WPP + 2*X*PLPP*QW-2*X 
+ *PLPP*W2 + X*QPP*ML2-X*QPP*W2 + X*WPP*ML2 
+ + X*WPP*W2 + PLPP*QW*ML2 + PLPP*QW*W2-2*QW*WPP 
+ *ML2-1./2.*QPP*ML4 + QPP*ML2*W2-1 . /2 . *QPP*W2**2) 

RI=RI + A2*C2*(-4*X**2*W2 + 4*X*QW*ML2 + 2*X* 
+ QW*W2-2*QW**2*ML2) 

SM=R1+RI/X+R2/X**2 

RETURN 

END 
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List of Tables 



Phenomenological values and source for the renormalized coupling constants| 



[M ). The quantities Y\ in the fourth column determine the scale depen-| 



dence ol" the ZTT/i] according to Eq. (16]). 



L r u and L\ 2 are not directly 



accessible to experiment. 



2 Occurrence of the low-energy coupling constants Li, . 


. , L\q and of the anomaly 


in the scmilcptonic decays discussed in chapter 2 





1.1 The quantities Xj,Nj. SD"^ and INT ± are evaluated with full phase space. 
IB with restricted kinematics (1.200 



1.2 Measured branching ratios T(K — > lui r y)/T t . 



The K e2 ry data are from ]5| , 



the data from 



3|. The last column corresponds ]3| to the cut QT3 

1.3 Chiral prediction at order for the branching ratios T(K — » lui'-f) /T t ' € 
cut used in the last column is given in Eq. (1.20 ) 



TEe) 



2.1 'l'heoretical values for the branching ratios for the decay K + — » jj, + ve + e for 



various cuts. 



2.2 Theoretical values for the branching ratios for the decay K~ 



e T re T e 



for 



various cuts. 



3.1 Kates of A; 3 decays. The number of events in the third column corresponds] 



to those data which are of relevance for the determination of the slope Aq ol 
the scalar form factor 



16 
20 
21 
27 

28 

38 

43 
43 



4.1 Branching ratios for tree level amplitudes for ~EZ > 30MeV and 0^ > 20° in 



the K rest frame 



4.2 Experimental results for decays 



4.3 



Coefficients for the loop amplitudes corresponding to the diagrams T 



1, 2, 3 in Fig. 4.2|. All coefficients c\ must be divided by 6\^2F 2 45 



4.4 Coefficients for the KjtZ loop amplitudes corresponding to the diagrams 7 
1, 2, 3 in Fig. 4.2 . All coefficients c\ must be divided by 6\^2F 2 



4.5 Branching ratios and expected number of events at DAFNE for K 



12a. 



4.6 Branching ratios and expected number of events at DAFNE for K\ 



46 
47 
47 
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5.1 Rates of K e4 decays 



5.2 Predictions of chiral symmetry following from the fit to the K e4 data 


441 


alone (column 3) and the combined determination from tttt 


55 


| and K e ± data 



(last column). The first column gives the prediction of the leading order 
term in the low-energy expansion of the tttt amplitude 
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List of Figures 



1.1 Contour plots lor / IB ,. . . , / INT ± 




. The numbering on the lines points] 


towards increasing modulus. The normalization is arbitrary 



1.2 Contour plots lor f m ,. . . , f mT ± 




. The numbering on the lines points] 


towards increasing modulus. The normalization is arbitrary 



1.3 The rate dP(x)/dx in (1.39] ) , evaluated with the form factors ( 1.38 ) and N tot = 
9 • 10 9 . The solid line corresponds to M K * = 890 MeV, M Kl = 1.3 GeV. The 
dashed line is evaluated with Mk* = 890 MeV, Mk x = oo and the dotted 
line corresponds to Mk* = Mk x = oo. The total number of events is also 
indicated in each case 

3.1 The vector and scalar form factors /+(£) and fp(t). 



3.2 The normalized slopes of the vector and the scalar form factors. Curve 1: the] 



normalized slope M%+df+(t)/dt. Curve 2: the normalized slope M%+dfo(t)/dt. 


Near the irK threshold to 


= {Mk + M^Yi the vector form factor 


behaves as] 


u(t) = u(t n ) + o 


(t-t ) 


, whereas f (t) = fo(t ) + 


Wt-h) 


. The slope] 



of the scalar form factor is therefore singular at t = {Mk + M^ 



4.1 Diagrammatic representation of the amplitude 



4.2 hoop diagrams (without tadpoles) for K/3 at Q{p i ). For K;3 7 , the photon 



must be appended on all charged lines and on all vertices 



FTT Kinematic variables lor K14 decays. TEe angle 6^ is dehned in Sj^T^ in Y. 



an 



d <b in S 



K ■ 



[5.2 The form factors f s (s n ) and g(s 7T ) (Eq. 5.64f) according to the chiral rep- 
resentation (electron mode). The dotted lines show the lowest order result 
(|5.47| ), and the dashed lines correspond to Li,L 2 and L 3 from (|5.63j) . The 
experimental result ( |5.37| ) is displayed by a solid line 
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